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ABSTRACT 
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studied  . 


TABLE  OF  CONTENTS 


I.  INTRODUCTION 15 

II.  EFFECT  OF  HEAT  INPUT  ON  COMBUSTION 19 

A.  INTRODUCTION  AND  DESCRIPTION  19 

B.  PROCEDURE 21 

C.  RESULTS 22 

1.  CASE  II-l   SQ  =  18000  BTU/ft.sq  hr  ....  22 

2.  CASE  II-2   SQ  =  20000  BTU/ft.sq  hr  ....  23 

3.  CASE  II-3   SQ  =  30000  BTU/ft.sq  hr  ....  28 

4.  CASE  II-4   SQ  =  40000  BTU/ft.sq  hr  ....  34 

D.  SUMMARY 40 

1.  Power  Relation 40 

2.  Relation  of  Temperature  and  Oxygen 
Concentration   40 

3.  Equilibrium  Temperature  for  Combustion  .  .  46 

III.  EFFECT  OF  THICKNESS  ON  COMBUSTION 49 

A.  INTRODUCTION 49 

B.  PROCEDURE 50 

C.  RESULTS 50 

1.  CASE  III-l   Thickness  =0.25  inches   ...  51 

2.  CASE  III-2   Thickness  =0.50  inches   ...  51 

3.  CASE  III-3   Thickness  =0.75  inches   ...  59 

4.  CASE  III-4   Thickness  =  1.00  inches   ...  61 

5.  CASE  III-5   Thickness  =2.00  inches   ...  64 

6.  CASE  III-6   Thickness  =  4.00  inches   ...  70 

7.  CASE  III-7   Thickness  =  6.00  inches   ...  73 

D.  SUMMARY 81 

1.  Power  Relation 81 

2.  Combustion  Start  Time 83 


IV.  PERMEABILITY 8  9 

A.  INTRODUCTION 8  9 

B.  PROCEDURE 91 

C.  RESULTS 92 

1.  CASE  IV- 1  Permeability  =  0.00015  ft2 
(d=D=0.0025  in) 92 

2.  CASE  IV- 2  Permeability  =  0.00058  ft2 
(d=D=0.0050  in) 93 

3.  CASE  IV- 3  Permeability  =  0.00233  ft2 
(d=D=0.0100  in) 97 

4.  CASE  IV-4  Permeability  =  0.00933  ft2 
(d=D=0.0200  in) 100 

D.  SUMMARY 103 

1.  Power  Relation 103 

2.  Combustion  Speed  105 

V.  REACTION  ORDER Ill 

A.  INTRODUCTION Ill 

B.  PROCEDURE Ill 

C.  RESULTS 112 

1.  CASE  V-l   Reaction  Order  =  0.25   ....  112 

2.  CASE  V-2   Reaction  Order  =  0.50   ....  114 

3.  CASE  V-3   Reaction  Order  =  0.75   ....  117 

4.  CASE  V-4   Reaction  Order  =1.00   ....  121 

D.  SUMMARY 123 

1.   Power  Relation 123 

VI.  CONCLUSIONS 132 

LIST  OF  REFERENCES 134 

INITIAL  DISTRIBUTION  LIST  135 


LIST  OF  TABLES 


I  VARIED  PARAMETERS  DURING  THE  PROCESS  FOR  SQ  = 

18000  BTU/FT.SQ  HR    27 

II  VARIABLE  PARAMETERS  DURING  THE  PROCESS  FOR  SQ 

=  20000  BTU/FT.SQ  HR    29 

III  VARIED  PARAMETERS  DURING  THE  PROCESS  FOR  SQ  = 

30000  BTU/FT.SQ  HR 35 

IV  VARIED  PARAMETERS  DURING  THE  PROCESS  FOR  SQ  = 

40000  BTU/FT.SQ  HR    41 

V  VARIED  PARAMETERS  DURING  THE  PROCESS  FOR 
THICKNESS  =0.25  INCHES   55 

VI  VARIED  PARAMETERS  DURING  THE  PROCESS  FOR 
THICKNESS  =0.50  INCHES   59 

VII  VARIED  PARAMETERS  DURING  THE  PROCESS  FOR 
THICKNESS  =0.75  INCHES   64 

VIII  VARIED  PARAMETERS  DURING  THE  PROCESS  FOR 
THICKNESS  =  1.00  INCHES 68 

IX  VARIED  PARAMETERS  DURING  THE  PROCESS  FOR 
THICKNESS  =2.00  INCHES   73 

X  VARIED  PARAMETERS  DURING  THE  PROCESS  FOR 
THICKNESS  =4.00  INCHES   77 

XI  VARIED  PARAMETERS  DURING  THE  PROCESS  FOR 
THICKNESS  =6.00  INCHES   82 

XII  RESULTS  FOR  ALL  THICKNESS  CASES 88 

XIII  RESULTS  FOR  PERMEABILITY    104 

XIV  POWER  RELATION  RESULTS    108 

XV  COMBUSTION  AND  EXTINCTION  TEMPERATURE  FOR 

REACTION  ORDER  CASES  129 


LIST  OF  FIGURES 


2.1  Temperature  vs  X/L  and  time  for  SQ  =  18000 
BTU/ft.sq  hr  and  TQ  =  46  seconds 23 

2.2  Temperature  vs  X/L  and  time  for  SQ  =  18000 
BTU/ft.sq  hr  and  TQ   =  45  seconds 23 

2.3  Oxygen  concentration  vs  X/L  and  time  for  SQ  = 

18000  BTU/ft.sq  hr  and  TQ  =  46  seconds 24 

2.4  Oxygen  concentration  vs  X/L  and  time  for  SQ  = 

18000  BTU/ft.sq  hr  and  TQ  =  45  seconds 24 

2.5  Reaction  rate  vs  X/L  and  time  for  SQ  =18000 
BTU/ft.sq  hr  and  TQ  =■  46  seconds 25 

2.6  Reaction  rate  vs  X/L  and  time  for  SQ  =  18000 
BTU/ft.sq  hr  and  TQ  =  45  seconds 25 

2.7  Combustion  and  extinction  carbon  temperature 

for  position  X/L  =  0 26 

2.8  Temperature  vs  X/L  and  time  for  SQ  =  20000 
BTU/ft.sq  hr  and  TQ  =  35  seconds 30 

2.9  Temperature  vs  X/L  and  time  for  SQ  =  20000 
BTU/ft.sq  hr  and  TQ  =  34  seconds 30 

2.10  Oxygen  concentration  vs  X/L  and  time  for  SQ  = 

20000  BTU/ft.sq  hr  and  TQ  =  35  seconds 31 

2.11  Oxygen  concentration  vs  X/L  and  time  for  SQ  = 

20000  BTU/ft.sq  hr  and  TQ  =  34  seconds 31 

2.12  Reaction  rate  vs  X/L  and  time  for  SQ  =  20000 
BTU/ft.sq  hr  and  TQ  =  35  seconds 32 

2.13  Reaction  rate  vs  X/L  and  time  for  SQ  =20000 
BTU/ft.sq  hr  and  TQ  =  34  seconds 32 

2.14  Extinction  and  combustion  carbon  temperature 

for  position  X/L  =  0 33 

2.15  Temperature  vs  X/L  and  time  for  SQ  =  30000 
BTU/ft.sq  hr  and  TQ  =  15  seconds 36 

2.16  Temperature  vs  X/L  and  time  for  SQ  =  30000 
BTU/ft.sq  hr  and  TQ  =  14  seconds 36 

2.17  Oxygen  concentration  vs  X/L  and  time  for  SQ 

=  30000  BTU/ft.sq  hr  and  TQ  =  15  seconds 37 

2.18  Oxygen  concentration  vs  X/L  and  time  for  SQ  = 

30000  BTU/ft.sq  hr,and  TQ  =  14  seconds  ......    37 

2.19  Reaction  rate  vs  X/L  and  time  for  SQ  =  30000 
BTU/ft.sq  hr  and  TQ  =  15  seconds 38 


2.20  Reaction  rate  vs  X/L  and  time  for  SQ  =  30000 
BTU/ft.sq  hr  and  TQ  =  14  seconds 38 

2.21  Extinction  and  combustion  carbon  temperature 

for  position  X/L  =  0 39 

2.22  Temperature  vs  X/L  and  time  for  SQ  =  40000 
BTU/ft.sq  hr  and  TQ  =  9  seconds 42 

2.23  Temperature  vs  X/L  and  time  for  SQ  =  40000 
BTU/ft.sq  hr  and  TQ  =  8  seconds 42 

2.24  Oxygen  concentration  vs  X/L  and  time  for  SQ  = 

40000  BTU/ft.sq  hr  and  TQ  =  9  seconds 43 

2.25  Oxygen  concentration  vs  X/L  and  time  for  SQ  = 

40000  BTU/ft.sq  hr  and  TQ  =  8  seconds 43 

2.26  Reaction  rate  vs  X/L  and  time  for  SQ  =  40000 
BTU/ft.sq  hr  and  TQ  =  9  seconds 44 

2.27  Reaction  rate  vs  X/L  and  time  for  SQ  =  40000 
BTU/ft.sq  hr  and  TQ  =  8  seconds 44 

2.28  Extinction  and  combustion  carbon  temperature 

for  position  X/L  =1 45 

2.29  Heat  flux  versus  combustion  temperature  (TQ  ) 
rectangular  plot 47 

2.30  Heat  flux  (SQ)  vs  combustion  temperature (TQ  ) 

Log  log  plot ...  48 

3.1  Temperature  vs  X/L  and  time  for  thickness  = 

0.25  in  Initial  carbon  temperature  =  1100  F  ....  52 

3.2  Temperature  vs  X/L  and  time  for  thickness  = 

0.25  in  Initial  carbon  temperature  =  1090  F  ....  52 

3.3  Oxygen  concentration  vs  X/L  and  time  for 
thickness  =  0.25  in  Initial  carbon  temperature 

=  1100  F 53 

3.4  Oxygen  concentration  vs  X/L  and  time  for 
thickness  =  0.25  in  Initial  carbon  temperature 

=  1090  F 53 

3.5  Reaction  rate  vs  X/L  and  time  for  thickness  = 

0.25  in  Initial  carbon  temperature  =  1100  F  ....  54 

3.6  Reaction  rate  vs  X/L  and  time  for  thickness  = 

0.25  in  Initial  carbon  temperature  =  1090  F  ....  54 

3.7  Temperature  vs  X/L  and  time  for  thickness  = 

0.50  in  Initial  carbon  temperature  =  960  F   ....  56 

3.8  Temperature  vs  X/L  and  time  for  thickness  = 

0.50  in  Initial  carbon  temperature  =  950  F   ....  56 

3.9  Oxygen  concentration  vs  X/L  and  time  for 
thickness  =  0.50  in  Initial  carbon  temperature 

=  960  F . 57 

3.10  Oxygen  concentration  vs  X/L  and  time  for 
thickness  =  0.50  in  Initial  carbon  temperature 

=  950  F 57 


3.11 
3.12 
3.13 
3.14 
3.15 

3.16 

3.17 
3.18 
3.19 
3.20 
3.21 

3.22 

3.23 
3.24 
3.25 
3.26 
3.27 

3.28 

3.29 
3.30 


Reaction  rate  vs  X/L  and  time  for  thickness  = 

0.50  in  Initial  carbon  temperature  =  960  F   ....  58 

Reaction  rate  vs  X/L  and  time  for  thickness  = 

0.50  in  Initial  carbon  temperature  =  950  F   ....  58 

Temperature  vs  X/L  and  time  for  thickness  = 

0.75  in  Initial  carbon  temperature  =  900  F   ....  60 

Temperature  vs  X/L  and  time  for  thickness  = 

0.75  in  Initial  carbon  temperature  =  890  F   ....  61 

Oxygen  concentration  vs  X/L  and  time  for 
thickness  =  0.75  in  Initial  carbon  temperature 
=  900  F 

Oxygen  concentration  vs  X/L  and  time  for 
thickness  =  0.75  in  Initial  carbon  temperature 
=  890  F 


Reaction  rate  vs  X/L  and  time  for  thickness  = 
0.75  in  Initial  carbon  temperature  =  900  F 

Reaction  rate  vs  X/L  and  time  for  thickness  = 
0.75  in  Initial  carbon  temperature  =  890  F 

Temperature  vs  X/L  and  time  for  thickness  = 
1.00  in  Initial  carbon  temperature  =  850  F 

Temperature  vs  X/L  and  time  for  thickness  = 
1.00  in  Initial  carbon  temperature  =  840  F 

Oxygen  concentration  vs  X/L  and  time  for 
thickness  =  1.00  in  Initial  carbon  temperature 
=  850  F .     .     . 

Oxygen  concentration  vs  X/L  and  time  for 
thickness  =  1.00  in  Initial  carbon  temperature 
=  840  F .  .  . 


Reaction  rate  vs  X/L  and  time  for  thickness  = 
1.00  in  Initial  carbon  temperature  =  850  F 

Reaction  rate  vs  X/L  and  time  for  thickness  = 
1.00  in  Initial  carbon  temperature  =  840  F 

Temperature  vs  X/L  and  time  for  thickness  = 
2.00  in  Initial  carbon  temperature  =  760  F 

Temperature  vs  X/L  and  time  for  thickness  = 
2.00  in  Initial  carbon  temperature  =  750  F 

Oxygen  concentration  vs  X/L  and  time  for 
thickness  =  2.00  in  Initial  carbon  temperature 
=  760  F .  .  . 

Oxygen  concentration  vs  X/L  and  time  for 
thickness  =  2.00  in  Initial  carbon  temperature 
=  750  F 

Reaction  rate  vs  X/L  and  time  for  thickness  = 
2.00  in  Initial  carbon  temperature  =  760  F 

Reaction  rate  vs  X/L  and  time  for  thickness  = 
2.00  in  Initial  carbon  temperature  =  750  F 


61 

62 
62 
63 
65 
65 

66 

66 
67 
67 
69 
70 

70 

71 
71 
72 


3.31  Temperature  vs  X/L  and  time  for  thickness  = 

4. 00  in  Initial  carbon  temperature  =  690  F   ....  74 

3.32  Temperature  vs  X/L  and  time  for  thickness  = 

4.00  in  Initial  carbon  temperature  =  680  F   ....  74 

3.33  Oxygen  concentration  vs  X/L  and  time  for 
thickness  =  4.00  in  Initial  carbon  temperature 

=  690  F 75 

3.34  Oxygen  concentration  vs  X/L  and  time  for 
thxckness  =  4.00  in  Initial  carbon  temperature 

=  680  F 75 

3.35  Reaction  rate  vs  X/L  and  time  for  thickness  = 

4.00  in  Initial  carbon  temperature  =  690  F   ....  76 

3.36  Reaction  rate  vs  X/L  and  time  for  thickness  = 

4.00  in  Initial  carbon  temperature  =  680  F   ....  76 

3.37  Temperature  vs  X/L  and  time  for  thickness  = 

6.00  in  Initial  carbon  temperature  =  650  F   ....  78 

3.38  Temperature  vs  X/L  and  time  for  thickness  = 

6.00  in  Initial  carbon  temperature  =  640  F   ....  79 

3.39  Oxygen  concentration  vs  X/L  and  time  for 
thickness  =  6.00  in  Initial  carbon  temperature 

=  650  F 79 

3.40  Oxygen  concentration  vs  X/L  and  time  for 
thickness  =  6.00  in  Initial  carbon  temperature 

=  640  F 80 

3.41  Reaction  rate  vs  X/L  and  time  for  thickness  = 

6.00  in  Initial  carbon  temperature  =  650  F   ....  80 

3.42  Reaction  rate  vs  X/L  and  time  for  thickness  = 

6.00  in  Initial  carbon  temperature  =  640  F   ....  81 

3.43  Thickness  versus  combustion  temperature 

Retangular  plot 84 

3.44  Thickness  versus  combustion  temperature  Log- log 

plot 85 

3.45  Thickness  vs  time  for  oxygen  concentration  at 

X/L  =  1.  to  reach  zero  Retangular  plot 86 

3.46  Thickness  vs  time  for  oxygen  concentration 

reach  zero  log-log  plot 87 

4.1  Geometry  of  a  typical  cell 90 

4.2  Temperature  vs  X/L  and  time  for  permeability  = 
0.00015  ft*  initial  carbon  temperatiure  =  750  F  .  .  93 


4.3 


4.5 


Temperature  vs  X/L  and  time  for  permeability  = 
0.00015  ftz  initial  carbon  temperatiure  =  750  F  .  .  93 


4.4    Oxygen  concentration  vs  XAL   and  time  for 
permeability  =  0.00015  ftz  initial 


y  =  0.00015  ft"1  initial  carbon 
temperatiure  =  760  F 94 

Oxygen  comcentration  vs  XAL    and  time  for 
permeability  =  0.00015  ft   initial  carbon 
temperatiure  =  750  F , 94 


10 


4.9 
4.10 


4.6  Reaction  rate  vs  X/L  and  time  for  permeability 

=  0.00015  ft   initial  carbon  temperatiure  =  760 

F 95 

4.7  Reaction  rate  vs  X/L  and  time  for  permeability 

=  0.00015  ft   initial  carbon  temperatiure  =  750 

F 95 

4.8  Temperature  vs  X/L  and  time  for  permeability  = 
0.00058  ftz  initial  carbon  temperatiure  =  850  F  .  .  96 

Temperature  vs  X/L  and  time  for  permeability  = 
0.00058  ftz  initial  carbon  temperatiure  =  840  F  .  .  97 

Oxygen  concentration  vs  XXL  and  time  for 
permeability  =  0.00058  ft   initial  carbon 
temperatiure  =  850  F 97 

Oxygen  concentration  vs  XXL  and  time  for 
permeability  =  0.00058  ft   initial  carbon 
temperatiure  =  840  F 98 

Reaction  rate  vs  X/L  and  time  for  permeability 

=  0.00058  ft   initial  carbon  temperatiure  =  850 

F 98 

Reaction  rate  vs  X/L  and  time  for  permeability 

=  0.00058  ft   initial  carbon  temperatiure  =  840 

F 99 

Temperature  vs  X/L  and  time  for  permeability  = 
0.00233  ftz  initial  carbon  temperatiure  =  980  F  .   100 

4.15   Temperature  vs  X/L  and  time  for  permeability  = 

0.00233  ftz  initial  carbon  temperatiure  =  970  F  .   100 


4.11 
4.12 
4.13 
4.14 


4.16 
4.17 


Oxygen  concentration  vs  XXL  and  time  for 
permeability  =  0.00233  ft   initial  carbon 
temperatiure  =  980  F 101 

Oxygen  comcentration  vs  X/L  and  time  for 
permeability  =  0.00233  ft   initial  carbon 
temperatiure  =  970  F 101 

4.18  Reaction  rate  vs  X/L  and  time  for  permeability 

=  0.00233  ftz  initial  carbon  temperatiure  =  980 

F 102 

4.19  Reaction  rate  vs  X/L  and  time  for  permeability 

=  0.00233  ft   initial  carbon  temperatiure  =  970 

F 102 

4.20  Permeability  versus  combustion  temperature 
retangular  plot 106 

4.21  Permeability  versus  combustion  temperature  log 

log  plot .   107 

4.22  Permeability  versus  t*  Retangular  plot   110 

4.23  Permeability  versus  t*  Log-log  plot Ill 

5.1  Temperature  vs  X/L  and  time  for  reaction  order 

=  0.25  Initial  carbon  temperature  =  980  F  .  .  .  .   114 

5.2  Temperature  vs  X/L  and  time  for  reaction  order 

=  0.25  Initial  carbon  temperature  =  970  F*.  .  .  .   115 


11 


5.3  Oxygen  concentration  vs  X/L  and  time  for 
reaction  order  =  0.25  Initial  carbon 

temperature  =  980  F 115 

5.4  Oxygen  concentration  vs  X/L  and  time  for 
reaction  order  =  0.25  Initial  carbon 

temperature  =  970  F 116 

5.5  Reaction  rate  vs  X/L  and  time  for  reaction 

order  =  0.25  Initial  carbon  temperature  =  980  F  .   116 

5.6  Reaction  rate  vs  X/L  and  time  for  reaction 

order  =  0.25  Initial  carbon  temperature  =  970  F  .   117 

5.7  Temperature  vs  X/L  and  time  for  reaction  order 

=  0.50  Initial  carbon  temperature  =  850  F  .  .  .  .   118 

5.8  Temperature  vs  X/L  and  time  for  reaction  order 

=  0.50  Initial  carbon  temperature  =  840  F  .  .  .  .   118 

5.9  Oxygen  concentration  vs  X/L  and  time  for 
reaction  order  =  0.50  Initial  carbon 

temperature  =  850  F 119 

5.10  Oxygen  concentration  vs  X/L  and  time  for 
reaction  order  =  0.50  Initial  carbon 

temperature  =  840  F 119 

5.11  Reaction  rate  vs  X/L  and  time  for  reaction 

order  =  0.50  Initial  carbon  temperature  =  850  F  .   120 

5.12  Reaction  rate  vs  X/L  and  time  for  reaction 

order  =  0.50  Initial  carbon  temperature  =  840  F  .   120 

5.13  Temperature  vs  X/L  and  time  for  reaction  order 

=  0.75  Initial  carbon  temperature  =  750  F  .  .  .  .   121 

5.14  Temperature  vs  X/L  and  time  for  reaction  order 

=  0.75  Initial  carbon  temperature  =  740  F  .  .  .  .   122 

5.15  Oxygen  concentration  vs  X/L  and  time  for 
reaction  order  =  0.75  Initial  carbon 

temperature  =  750  F 122 

5.16  Oxygen  concentration  vs  X/L  and  time  for 
reaction  order  =  0.75  Initial  carbon 

temperature  =  740  F 123 

5.17  Reaction  rate  vs  X/L  and  time  for  reaction 

order  =  0.75  Initial  carbon  temperature  =  750  F  .   123 

5.18  Reaction  rate  vs  X/L  and  time  for  reaction 

order  =  0.75  Initial  carbon  temperature  =  740  F  .   124 

5.19  Temperature  vs  X/L  and  time  for  reaction  order 

=  1.00  Initial  carbon  temperature  =  670  F  .  .  .  .   125 

5.20  Temperature  vs  X/L  and  time  for  reaction  order 

=  1.00  Initial  carbon  temperature  =  660  F  .  .  .  .   126 

5.21  Oxygen  concentration  vs  X/L  and  time  for 
reaction  order  =  1.00  Initial  carbon 

temperature  =  670  F 126 

5.22  Oxygen  concentration  vs  X/L  and  time  for 
reaction  order  =  1.00  Initial  carbon 

temperature  =  660  F * 127 


12 


5.23  Reaction  rate  vs  X/L  and  time  for  reaction 

rate=  1.00  Initial  carbon  temperature  =  670  F  .  .   127 

5.24  Reaction  rate  vs  X/L  and  time  for  reaction 

order  =  1.00  Initial  carbon  temperature  =  660  F  .   128 

5.25  Reaction  order  versus  temperature  Retangular 

plot 130 

5.26  Reaction  order  versus  temperature  Log-log  plot   .   131 

5.27  Reaction  order  versus  temperature  Semi  log  plot  .   132 


13 


ACKNOWLEDGEMENTS 

I  wish  to  express  my  sincere  appreciation  to  Dr.  David 
Salinas  for  his  guidance,  assistance  and  encouragement 
during  the  pursuit  of  this  study.  To  my  wife  and  our  sons, 
for  the  encouragement,  patience  and  understanding,  I  am 
deeply  grateful. 


14 


I.  INTRODUCTION 

The  goal  of  this  thesis  is  to  study  the  combustion 
behavior  of  porous  graphite  .  The  mathematical  model  devel- 
oped by  VATIKIOTIS  (1)  is  a  transient,  one-dimensional  model 
of  a  porous  medium  consisting  of  spherical  particles  ,  or 
cylindrical  fibers  of  graphite  . 

The  model  covers  several  aspects  :  mass  transfer,  heat 
transfer,  combustion  ,  airflow  ,  and  temperature  dependency 
of  thermophysical  properties.  Geometric  parameters  of  the 
porous  medium  include  porosity  (p)  and  medium  thickness  (L) . 
The  porous  medium  has  interconnected  pores  which  permit  air 
flow  through  the  medium.  This  air  flow  was  modeled  by  the 
continuity  equation  and  Darcy ' s  Law  .  The  porous  medium  is 
characterized  by  a  number  of  parameters  including  porosity  , 
specific  internal  area  (Z),  tortuosity  (r)  and  permeability 
(m) .  These  parameters  are  defined  as  follows.  Given  a  unit 
square  cell  of  dimension  D,  and  particle  diameter  d,  we  have 

p  =  1  -  v    I    4  (d/D)2  (1.1) 


z  =  1/2  ,7  d2/D3  (1.2) 


=  1.4  (1.3) 


m  =    Sp3/z2  (1.4) 


Porosity,'  p,  is  defined  as  the  ratio  of  void  volume  per 
volume  of  unit  all  .   The  specific  internal  area,  Z,  is  the 
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ratio  of  internal  surface  area  to  bulk  volume  .  The 
tortuosity  is  defined  as  the  ratio  of  the  length  of  the 
flowpath  of  a  fluid  particle  to  the  straight  line  distance. 

The  particle's  size  decreases  as  the  carbon  is  consumed 
and  all  geometric  properties  which  depend  on  particle  diam- 
eter are  functions  of  time  and  position  .  This  model 
considers  that  the  carbon  matrix  remains  rigid  as  the 
particle  diameter  decreases  and  thus  porosity  increases  with 
combustion. 

The  magnitude  of  REYNOLDS  number  defines  the  motion  of 
the  fluid  which  can  be  molecular  ,  viscous  or  inertial.  In 
the  case  of  porous  media  ,  the  flow  is  dominated  by  viscous 
and  inertial  effects.  Darcy ' s  Law  for  the  range  of  REYNOLD'S 
number  where  viscous  effects  dominate  is  the  equation  that 
governs  the  fluid  flow  in  porous  media.  Two  velocities  which 
distinguish  the  flow  through  the  medium  are  filter  velocity 
Q,  and  pore  velocity  u.  Filter  velocity  is  given  by  : 

Q  =  -  m/M  (  <9P/3x  )  (1.5) 

where  u  is  the  viscosity  of  the  fluid  and  P  is  the  pressure. 
The  pore  velocity  u  is  given  by  : 

U  =  Q/P  (1.6) 

The  hypothesis  of  the  Dupuit-Forcheimer  assumption  is 
that  the  local  pore  velocity  is  greater  than  the  filter 
velocity.  The  actual  velocity  in  a  pore  is  a  function  of 
position  within  the  pore  .  The  Dupuit-Forcheimer  relation 
defines  an   average  velocity  in  a  pore. 

The  model  of  combustion  adopted  was  the  model  of  N.  N. 
Semenov  (2).  The  relation  of  reaction  rate  to  temperature 
and  oxygen  concentration,  and  the  interaction  of  heat  gener- 
ation and  heat  transfer  are  fundamental  to  this  model  .  If 
the  heat  transfer  dominates,  extinction  will  occur.  If  the 
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heat  generation  dominates  then  combustion  will  result.  The 
reaction  rate  equation  was  adopted  from  Arrhenius ' s  Law  for 
a  simple  reaction  where  the  rate  depends  on  the  concentra- 
tion of  the  reactants  and  not  on  the  products.  The  heat 
generated  by  an  exothermal  reaction  is  obtained  by  multi- 
plying the  reaction  rate  by  the  heat  of  combustion  . 

A  study  of  the  curve  of  heat  generation  versus  tempera- 
ture presented  by  Frank-Kamenet skii  distinguishes  two  phases 
during  combustion.  During  the  initial  kinetic  phase  ,  the 
rate  of  reaction  and  the  temperature  are  lowest.  During  This 
phase  ,  there  is  an  excess  supply  of  oxygen.  In  the  kinetic 
regime  the  reation  rate  increases  exponentially  with 
increasing  temperature.  The  kinetic  regime  is  followed  by 
the  diffusion  regime.  This  phase  is  distinguished  by  high 
temperatures  and  reaction  rates  ,  and  the  reaction  is 
limited  by  a  lack  of  oxygen. 

It  is  assumed  that  the  chemical  reaction  for  the  combus- 
tion process  produces  carbon  monoxide  and  carbon  dioxide. 
The  ratio  of  the  mass  rates  of  carbon  monoxide  to  carbon 
dioxide  depends  on  the  temperature.  Increasing  temperature 
results  in  an  increase  of  this  CO  to  COo  ratio. 

The  model  of  heat  transfer  includes  three  basic  mecha- 
nisms :  convection  ,  conduction  and  radiation  .  The  heat 
transfer  equation  used  by  VATIKIOTIS  (1)  also  includes  a 
heat  generation  term  due  to  combustion.  All  thermophysical 
properties  were  treated  as  temperatures  dependent.  These 
properties  include  viscosity,  conductivity  and  density. 
Energy  balance  equations  were  constructed  for  control 
volumes  of  the  porous  medium  and  the  air  flowing  through  the 
porous  medium  . 

A  third  field  equation  of  the  model  was  obtained  by  a 
mass  balance  of  oxygen  .  This  equation  includes  molecular 
diffusion,  convective  transport  and  a  sink  term  to  account 
for  the  depletion  of  oxygen  due  to  combustion. 
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These  equations  resulted  in  a  system  of  coupled,  nonli- 
near transient  field  equation  .  These  field  equations, 
together  with  boundary  and  initial  conditions  ,  define  the 
problem.  They  were  solved  numerically  by  using  the  Galerkin 
finite  element  method.  Details  of  the  formulation  are 
presented  in  reference  1  . 
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II.  EFFECT  OF  HEAT  INPUT  ON  COMBUSTION 

A.   INTRODUCTION  AND  DESCRIPTION 

In  this  section  the  effect  of  heat  input  on  the  combus- 
tion behavior  of  porous  medium  was  investigated.  The  porous 
medium  is  subjected  to  a  heat  source  on  its  surface  .  Two 
parameters  of  main  importance  are  considered:  The  magnitude 
of  the  heat  flux  of  the  source  to  the  porous  medium  (SQ), 
and  the  time  that  the  heat  flux  remains  active  (TQ) . 

To  do  this  analysis,  all  parameters  of  the  system  are 
fixed  except  the  time  period  of  heating  (TQ)  and  the  amount 
of  heat  flux  (SQ) . 

The  purpose  of  this  section  is  to  determine  the  relation 
between  SQ  and  TQ  which  will  result  in  combustion.  This  is 
achieved  as  follows: 

For  a  specific  value  of  heat  flux  SQ,  the  program  is  run 
with  different  values  of  TQ  until  combustion  occurs  at  (TQ) 
and  extinction  occurs  at   (TQ)    ,  where<5  is  a  small  time 
period  .  Here  we  took<5  seconds  . 

TQC   =  TQe  +  8  (2.1) 

Four  cases  of  heat  flux  were  studied  in  this  section: 

CASE  II. 1  SQ  =  18000  BTU/ft.sq  hr 

CASE  II. 2  SQ  =  20000  BTU/ft.sq  hr 

CASE  II. 3  SQ  =  30000  BTU/ft.sq  hr 

CASE  II. 4  SQ  =  40000  BTU/ft.sq  hr 

The  fixed  parameters  used  in  all  cases  are: 

Ambient   temperature  =80   F 
Ambient  temperature  =  2117  lb/ft.sq 
Tortuosty   =  1.400 
Filament  diameter   =  0.0004167  ft 
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Thickness  of  matrix  laminate  =  0.0004167  ft 

Thickness  of  porous  medium  =  0.02083  ft 

Gas  constant  for  air  =  53.34  lbf  ft/lbm  R 

Conductivity  of  filament  =  86  BTU/lbm  H.F. 

Specific  heat  of  filament  =  0.703  lbm/Cf 

Emissivity  of  filament  =  0.90 

Shape  factor  for  int.HF.XFER  coefficient  =  1.00 

Characteristc  length  of  medium  =  1  inch 

Heat  of  reaction   =  14090  BTU/Cf 

Reaction  order  =  0.50 

Stochiometric  ratio  (fuel/air)  =  0.375 

Reaction  coefficient  =  2065000  lbm/Cf  H 

Activation  energy  coefficient  =  28840  deg  R 

Pressure  differential  across  thickness  =  50  lb/ft. sq 

Initial  carbon  and  air  temperature  =80   F 

Initial  uniform  oxygen  concentration  =  0.0172  lbm/ft.cu 

The  equations  used  in  this  program  permit  three  sets  of 
boundary  conditions  that  are  approximations  to  physical 
situations.  The  boundary  conditions  selected  are  the  best 
approximations  of  a  one  dimensional  model.  The  set  of 
boundary  conditions  used  in  this  case  are  typical  of  thermal 
flow  reactor  without  radiation  from  the  boundary  surfaces  . 
for  the  problems  of  this  investigation,  the  following 
boundary  conditions  were  used  : 

-  No  heat  transfer  from  either  end  of  the  porous  medium, 
dt/dx  =0  at   X/L  =  0   and   X/L  =  1        (2.2) 


-  Danckwerts'  boundary  conditions  :  Conditions  for  the 
air  temperature  and  oxygen  concentration  were  assumed.  This 
means  : 
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Ka  dTa/3x  =  Pacau(Ta-Too)  X/L  =  0       (2.3) 


dTa/dx  =0  X/L  =  1       (2.4) 


De  d<£  /ax  =  u  (^  -  0<oq)  X/L  =  0       (2.5) 


(Ke  +  Kr)  aTc/ax  =  -a  e(Tc4  -  T4^)       X/L  =  1       (2.6) 


KadTa/dx  =  p.acau(Ta  -  Too)  X/L  =  1       (2.7) 


dTa/dx  -  0  X/L  =  1       (2.8) 


Ded0/dx  =  u(0  -  0^)  X/L  =  1       (2.9) 

B .   PROCEDURE 

For  each  case  ,  the  amount  of  heat  flux  (SQ)  was  fixed 
and  several  values  of  shut  off  time  (TQ)  were  chosen  until 
values  of  TQC  and  TQe  were  found  such  that  : 

TQ,  »  TQ,  (2.10) 


TQC  =  TQe  +  8  (2.11) 

where  5  is  a  small  number. 
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C.   RESULTS 

The  results  for  each  case  were  obtained  in  both  numer- 
ical and  graphical  form  .  The  behavior  of  carbon  tempera- 
ture, oxygen  concentration  and  reaction  rate  for  combustion 
and  extinction  are  shown  as  surfaces  versus  time  and  posi- 
tion The  program  also  calculates  the  air  temperature  in  the 
porous  medium.  However  the  air  temperature  surfaces  are  not 
shown  as  they  are  quite  similar  ,  in  most  cases  to  carbon 
temperature . 

There  are  three  surface  graphics  for  each  heat  flux 
case : 

-Carbon  temperature  versus  position  and  time 

-Oxygen  concentration  versus  position  and  time 

-Reaction  rate  versus  position  and  time 

In  adition  ,  there  are  two  curves  for  each  heat  flux 
(SQ)  ,  one  showing  the  variation  of  carbon  temperature 
versus  time  for  combustion  and  the  other  showing  the  varia- 
tion of  carbon  temperature  for  extinction  .  These  is  two 
curves  are  for  the  point  x/L=0.  (i.e.  where  the  air  enters 
the  porous  medium)  . 

1.   CASE  II-l   SQ  =    18000  BTU/ft.sq  hr 

The  values  for  shut  off  time  bounding  extinction  and 
combustion  are  : 

for  combustion   TQ   =  46  seconds 
for  extinction   TQ   =  47  seconds 

During  the  transient  analyses  ,  system  parameters 
varied  .  for  each  analysis  ,  the  maximum  and  minimum  values 
of  some  of  these  parameters  provide  some  insight  into  the 
character  of  the  particular  case  .  These  parameters  and 
their  values  are  given  in  Table  I 

The  graphical  results  are  shown  in  Figures  2.1  to 
2.6 
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Figure  2.1    Temperature  vs  X/L  and  time 
for  SQ  =  18000  BTU/ft.sq  hr  and  TQ  =  46  seconds 
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Figure  2.2    Temperature  vs  X/L  and  time 
for  SQ  =  18000  BTU/ft.sq  hr  and  TQ   =  45  seconds 
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OXYGEN  CON:.  5URFRCE  FRQ^l  GRRf3E 


Figure  2.3    Oxygen  concentration  vs  X/L  and  time 


for  SQ  =  18000  BTU/ft.sq  hr  and  TQ  =  46  seconds 
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TGEN  CONG,  sur-rcl  from  g=rf3: 


Figure  2.4    Oxygen  concentration  vs  X/L  and  time 
for  SQ  =  18000  BTU/ft.sq  hr  and  TQ  =  45  seconds. 
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PlBCTION  PPTl  5URFP.CE  rPCT 


Figure    2.5         Reaction    rate    vs    X/L    and    time 
for    SQ    =18000    BTU/ft.sq    hr    and    TQ    =    46    seconds 

p:p::]o;j  ppt:  su^p::  rPC^  G=Pr3: 


Figure  2.6    Reaction  rate  vs  X/L  and  time 
for  SQ  =  18000  BTU/ft.sq  hr  and  TQ  =  45  seconds 
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TIME  VERSUS  GRAPHITE  TEMPERATURE 


TIME  (S) 
FOR  NUMBER  POINT  1  AND  SQ=i8000 


Figure  2.7    Combustion  and  extinction  carbon  temperature 

for  position  X/L  =  0.. 
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2.   CASE  II-2   SQ  =  20000  BTU/ft.sg  hr 

For  this  case  the  heat  flux  duration  times  for  the 
conditions  of  extinction  and  combustion  are  : 


for  combustion     TQ   =  35  seconds 
for  extinction     TQ   =  34  seconds 


During  the  transient  history  of  this  case  several 
parameters  varied.  The  maximum  and  minimum  of  these  changing 
variables  are  in  the  Table  II  .  It  is  necessary  to  emphasize 
that  these  extrema  results  do  not  occur  at  the  same  time  or 
same  position  (X/L) .  They  are  maximum  and  minimum  values, 
obtained  at  different  times  and  position  during  the  tran- 
sient process.  These  values  show  the  characteristc  behavior 
of  this  case  . 

The  behavior  of  carbon  temperature  ,  oxygen  concen- 
tration and  reaction  rate  are  shown  in  Figures  2.8  to  2.13  . 
These  surfaces  are  an  important  characteristic  of  this  case 
because  they  show  the  evolution  of  the  combustion  process  or 
extinction  process  for  each  position  during  the  transient 
time  . 

3.   CASE  II-3   SQ  =  30000  BTU/ft.sg  hr 

The  values  of  heating  time  that  result  in  combustion 
and  extinction  for  this  value  of  heat  flux  are  : 

for  combustion    TQC   =  15  seconds 
for  extinction    TQe   =  14  seconds 

During  the  transient  analysis  ,  some  variation  of 
values  occured  to  several  parameters.  The  limits  of  varia- 
tion ,  maximum  and  minimum  values  ,  are  shown  in  Table  III  . 
The  process  of  combustion  and  extinction  have  particular 
characteristics  that  determine  the  variation  of  the  values 
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Figure    2.8         Temperature    vs    X/L    and    time 
for    SQ    =    20000    BTU/ft.sq    hr    and    TQ    =    35    seconds 
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Figure    2.9         Temperature    vs    X/L    and    time 
for    SQ    =    20000    BTU/ft.sq    hr    and    TQ    =    34    seconds 
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OXYGEN  CONC.  SUPFRCE  rRO"  GRRF3E 


Figure  2.10    Oxygen  concentration  vs  X/L  and  time 
for  SQ  =  20000  BTU/ft.sq  hr  and  TQ  =  35  seconds. 

OXYGEN  CONC.  SURFACE  rR0n  GRRF3E 


Figure  2.11    Oxygen  concentration  vs  X/L  and  time 
for  SQ  =  20000  BTU/ft.sq  hr  and  TQ  =  34  seconds. 
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REACTION  RRi:  5URFRCE  rRQM  GR?F3E 


Figure  2.12    Reaction  rate  vs  X/L  and  time 
for  SQ  =  20000  BTU/ft.sq  hr  and  TQ  =  35  seconds 
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Figure  2.13    Reaction  rate  vs  X/L  and  time 
for  SQ  =20000  BTU/ft.sq  hr  and  TQ  =  34  seconds 
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TIME  VERSUS  GRAPHITE  TEMPERATURE 


40 

TIME  (S) 
FOR  NUMBER  POINT  1  AND  SQ=20000 


Fieure  2.14    Extinction  and  combustion  carbon  temperature 


for  position  X/L  =  0.. 
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of  these  parameters.  The  maximum  and  minimum  values  give  a 
general  idea  of  system  behavior  .  These  limiting  values  do 
not  occur  at  same  position  (X/L)  or  same  time  . 

A  better  idea  of  the  development  of  temperature  , 
oxygen  concentration  and  reaction  rate  can  be  obtained  from 
Figures  2.15  to  2.20  These  surfaces  show  the  values  of  these 
parameters  for  each  point  (X/L)  during  the  transient  time  . 

Figure  2.14  shows  that  at  X/L  =  0.  the  transient 
behavior  for  extinction  and  combustion  are  the  same  until 
TQe  =  14  seconds  ,  the  time  the  heat  flux  is  terminated  . 
After  this  time  the  temperature  for  the  extinction  case 
decreases  to  ambient  temperature  ;  while  for  the  combustion 
case  ,  the  temperature  increases  reaching  a  peak  temperature 
of  3200  F  at  28  seconds  and  thereafter  decreases  to  an  equi- 
librium combustion  temperature  of  about  2500  F  after  about 
50  seconds  (see  Figure  2.14). 


4.   CASE  II-4   SQ  =  40000  BTU/ft.sq  hr 

The  heat  flux  of  this  case  produced  combustion  and 
extinction  of  the  porous  medium  at  the  following  times  : 


for  combustion    TQC  =  9  seconds 
for  extinction    TQ   =  8  seconds 

During  the  transient  process  ,  the  output  of  the 
program  showed  parameters  th  at  didn't  change  with  time  and 
others  that  did  .  The  variation  of  these  parameter's  values 
show  the  development  of  the  combustion  and  extinction 
phenomenon  for  this  value  of  SQ  .  The  maximum  and  minimum 
values  of  several  parameters  are  shown  in  Table  IV  .  Again 
these  values  do  not  correspond  to  the  same  time  . 

The  characteristics  of  this  case  ,  carbon  tempera- 
ture ,  oxygen  concentration  and  reaction  rate  ,  as  function 
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TEMPERATURE  SURFACE  FROM  GRRF3E 


Figure    2.15         Temperature   vs    X/L   and    time 
for    SQ    =    30000    BTU/ft.sq    hr    and    TQ    =    15    seconds 

TEMPERATURE  SURfRCE  fRCM  GRRF3E 


Figure    2.16         Temperature    vs    X/L    and    time 
for    SQ    =    30000    BTU/ft.sq   hr    and    TQ    =    14    seconds 
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OXYGEN  CONC.  SURFACE  FROM  GRRF3E 


Figure  2.17    Oxigen  concentration  vs  X/L  and  time 
for  SQ  =30000  BTU/ft.sq  hr  and  TQ  =  15  seconds. 
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Figure  2.18    Oxigen  concentration  vs  X/L  and  time 
for  SQ  =  30000  BTU/ft.sq  hr , and  TQ  =  14  seconds. 
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REACTION  RATE  SURFACE  FROM  GRAF3E 


Figure  2.19    Reaction  rate  vs  X/L  and  time 
for  SQ  =  30000  BTU/ft.sq  hr  and  TQ  =  15  seconds 
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Figure  2.20    Reaction  rate  vs  X/L  and  time 
for  SQ  =  30000  BTU/ft.sq  hr  and  TQ  =  14  seconds 
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TIME  VERSUS  GRAPHITE  TEMPERATURE 
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Figure  2.21 


Extinction  and  combustion  carbon  temperature 
for  position  X/L  =  0.. 
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of  time  and  position  (X/L)  are  shown  as  graphic  surfaces  in 
Figures  2.22  to  2.27  .  These  graphical  results  give  an  idea 
about  the  behavior  of  each  point  of  the  transient  history  . 

For  X/L  =  0  ,  both  the  combustion  and  extinction 
cases  show  similar  behavior  to  this  time  equal  to  about  8 
seconds  ,  when  the  heat  flux  is  terminated  .  After  this 
time,  for  the  extinction  case  ,  the  temperature  decreases  to 
ambient  temperature  .  For  the  combustion  case  ,  the  temper- 
ature at  X/L  =  0  increased  to  about  3300  F  ,  and  thereafter 
decreased  to  around  2500  F  (temperature  of  combustion) 
Both  the  extinction  and  combustion  cases  achieve  new  equi- 
librium temperatures  ,  about  2500  F  for  the  combustion  cases 
and  ambient  temperature  for  the  extinction  cases  ,  after  35 
seconds  (see  Figures  2.28  ). 

D .   SUMMARY 

Here  ,  some  observations  are  made  about  the  effects  of 
heat  flux  input  and  duration  of  heat  flux  on  the  combustion 
process . 

1 .  Power  Relation 

For  each  case  ,  a  pair  of  values  (SQ,(TQ)  )  was 
obtained  .  A  plot  of  these  point  is  shown  in  Figures  2.29 
and  2.30  .  It  is  observed  that  the  relation  between  SQ  and 
TQ  of  these  graphics  yields  the  approximate  power  relation 
(equation  2 . 12 )  . 

2 .  Relation  of  Temperature  and  Oxygen  Concentration 
During  the  initial  heating  of  the  porous  medium  , 

the  temperature  increases  while  the  oxygen  is  being 
defleted.  After  the  heat  flux  is  removed  ,  either  extinction 
or  combustion  will  proceed.  If  extinction  occurs  ,  the 
carbon  temperature  decreases  to  ambient  temperature  while 
the  oxygen  is  being  restored  to  ambient  concentreit ion  . 

If  combustion  occurs  , then  the  oxygen  inside  the 
medium  goes  to  zero  and  remains  there   .   The  only  oxygen 


40 


CO 

H 

u. 

ZD 
H 
CQ 

O 
O 
O 
O 


co 

3 

E 
O 

a 


O 
m 


o     o 

O     r-> 


o      rn 
00     o 

H       O 


ON 


CN 


CN 


ON 

r»» 

o 

r>H 

r^ 

r«N 

o 

• 

ON 

<r 

• 

m 

m 

• 

CO 

CN 

ON 

>— i 

wo 

co 


o 

CO 

co 

u 
o 


3 

£ 

3      -H 

0       3 


U 

C 

•H  = 


OJ       — i 

c 


o 


O 
CO 


VO 
O 


CN 


O 


o 

CN 


ON 

cNi     m 

ON      m 


CN 


CN       CN 


ON 

--* 

O 

m 

oo 

O 

r—i 

<r 

^-1 

CN 

ON 

ao 

00 

ao 

-- 

m 

3 
U 


t 

o 


<4-l 


co 

t-l 

cu 
■u 
CU 
E 
CO 
U 

ca 
a 


c 

10 


J3 
ccj 
cu 
E 

at 
a 


c 
to 


U-|        ±j 


4-1 


0) 

u 

3 
CO 
V) 

cu 

u 

a. 


c 
cu 

•H 

T3 
C3 
U 

00 

HI 

l-l 

3 

01 

to 

01 

u 
a. 


u 
o 

I— I 

cu 

> 

cu 

u 
o 
a. 


u 
cu 

_a 
E 
3 
C 

co 

T3 

o 
c 
>> 
cu 
a5 


u 

0) 

> 

3 
O 

u 


41 


TEMPERATURE  SURFACE  FROM  GRRF3E 


Figure  2.22    Temperature  vs  X/L  and  time 
for  SQ  =  40000  BTU/ft.sq  hr  and  TQ  =  9  seconds 
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Figure  2.23    Temperature  vs  X/L  and  time 
for  SQ  =  40000  BTU/ft.sq  hr  and  TQ  -    8  seconds 
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OXYGEN  CONC.  SURFRCE  FROM  GRPF3E 
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Figure  2.24    Oxygen  concentration  vs  X/L  and  time 
for  SQ  =  40000  BTU/ft.sq  hr  and  TQ  =  9  seconds. 
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Figure  2.25    Oxigen  concentration  vs  X/L  and  time 
for  SQ  =  40000  BTU/ft.sq  hr  and  TQ  =  8  seconds. 
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RERCTION  RRTE  SURFRCl  FROM  GRRF3E 
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Figure  2.26    Reaction  rate  vs  X/L  and  time 
for  SQ  =  40000  BTU/ft.sq  hr  and  TQ  =  9  seconds 

RERCTION  RRTE  SURFRCE  FROn  SRRF3E 


Figure  2.27    Reaction  rate  vs  X/L  and  time 
for  SQ  =  40000  BTU/ft.sq  hr  and  TQ  =  8  seconds 
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TIME  VERSUS  GRAPHITE  TEMPERATURE 


•        < i      •  < »•   •   «-■ 


O*-  -©» 


< i-    •    i 1    ■■   >  • 


•<  •  •    ; t- 

;.  .    i — ;. 


:  i  <       --i       ••< i «■• 


■*■■ !--     ••<■ «• 

•■ »••     -■{■ i- 


LEGEND 
o    WITH  C0MBUST010N 


i.   .       ..  ..      ,  ...     ... |......«.. 

-e — e «       o  •   o     e  "■ 


40 


— T" 

to 


T"" 


70  M 


TIME  (S) 
FOR  NUMBER  POINT  1  AND  SQ=40000 


Figure  2.28 


Extinction  and  combustion  carbon  temperature 
for  position  X/L  =  1. 
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TQC  =  3.4  1010/(SQ)2.09  (2.12) 

present  is  at  the  entrance  at  X/L  =  0  .  The  temperature 
continues  to  increase  to  ' a  maximum  ,  and  thereafter  the 
temperature  decreases  .  Eventually  the  entire  medium  would 
achieve  a  uniform  combustion  temperature  . 

3 .   Equilibrium  Temperature  for  Combustion 

In  all  combustion  cases  ,  the  behavior  was  similar  . 
During  the  initial  combustion  period  ,  the  temperature  of 
the  medium  at  X/L  =  0  rose  between  2000  F  (CASE  II-l)  to 
3400F  (CASE  II-4)  .  At  this  time  ,  on  the  other  side  of  the 
medium  at  X/L  =1.  ,  the  temperature  of  the  medium  was 
approximately  1500  F  .  This  during  thus  early  transient 
stage  there  is  aAT  of  1300  F  (for  CASE  II-l)  and  1900  F 
(for  CASE  II-4)  across  the  medium  .  As  time  proceeds  the 
temperature  at  X/L  =  0  decreased  ,  while  at  X/L  =1.  ,  the 
temperature  increased  .  That  is  ,  with  progressing  time  , 
the  medium  moves  towards  a  uniform  temperature  .  In  the 
present  case  ,  this  equilibrium  temperature  of  combustion  is 
around  2400  F  to  2500  F  for  all  cases  . 

The  mechanism  of  equilibrium  of  the  temperature 
during  combustion  is  due  to  the  heat  transfer  by  air  .  As 
the  air  flows  through  the  porous  medium  it  transports  the 
heat  from  the  hot  front  part  of  the  medium  to  the  cooler 
back  part  ,  until  all  points  of  the  medium  have  the  same 
temperature  . 


46 


SQ  VERSUS   TQ 
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Figure  2.29    Heat  flux  versus  combustion  temperature  (TQ  ) 

rectangular  plot. 
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Figure  2.30 


Heat  flux  (SQ)  vs  combustion  temperature (TQC ) 
Log  log  plot . 


48 


III.  EFFECT  OF  THICKNESS  ON  COMBUSTION 

A.   INTRODUCTION 

The  porous  medium  has  several  geometric  parameters  that 
can  have  direct  influence  on  the  combustion  process.  The 
effect  of  the  thickness  parameter  on  system  behavior  is 
analysed  in  this  section.  All  other  parameters  are  fixed 
except  the  thickness  of  the  medium. 

The  purpose  of  this  section  is  to  determine  the  relation 
of  the  thickness  on  the  minimum  temperature  to  start  combus- 
tion (i.e.,  the  ignition  temperature). 

For  a  specified  value  of  thickness,  the  program  is  run 
with  different  values  of  initial  temperature.  This  initial 
temperature  is  equal  at  all  points  of  the  medium.  The  lowest 
value  of  initial  temperature  that  results  in  combustion  is 
the  ignition  temperature  Tc.  Extinction  occurs  at  Tc-<5  , 
wherefi  is  a  small  increment  of  temperature.  The  value  of<5 
used  in  this  parametric  analysis  was  10  F  . 

Seven  cases  with  different  thicknesses  were  studied  in 
this  section: 

CASE  III-l  THICKNESS  =0,25  inches 
CASE  III-2  THICKNESS  =0,50  inches 
CASE  III-3  THICKNESS  =0,75  inches 
CASE  III-4  THICKNESS  =  1.00  inches 
CASE  III-5  THICKNESS  =2.00  inches 
CASE  III- 6  THICKNESS  =4.00  inches 
CASE  III- 7  THICKNESS  =6.00  inches 
Each  case  is  characterised  by  its  thickness.  The  fixed 
parameters  used  in  all  cases  are: 

-  Ambient  temperature  =  80  F 

-  Ambient  pressure  =  2117  lb/ft.  sq. 

-  Tortuosity  =  1.400 

-  Filament  diameter  =  0.0004167  ft. 
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-  Thickness  of  matrix  laminate  =  0.0004167  ft. 

-  Thickness  of  porous  medium  =  0.02083  ft. 

-  Gas  constant  for  air  =  53.34  lbf.ft/lbm.R 

-  Conductivity  of  filament  =  86  BTU/lbf  HF 

-  Specific  heat  of  filament  =  0.703  lbm/cf. 

-  Emissivity  of  filament  =  0.90 

-  Shape  factor  for  int.HF.XFER  coefficient  =  1.00 

-  Heat  of  reaction  =  14090  BTU/CF 

-  Reaction  order  =  0.50 

-  Stoichiometric  ratio  (full/air)  =  0.375 

-  Reaction  rate  coefficient  =  2065000  lbm/CF  hr 

-  Activation  energy  coefficient  =  28.840  deg  R 

-  Pressure  differential  across  thickness  =  50  lb/ft 

-  Initial  uniform  oxygen  concentration  =  0.0172  lbm/ft 

B .  PROCEDURE 

For  each  case,  the  characteristic  length  of  the  medium 
(thickness)  was  fixed  and  values  of  the  uniform  initial 
temperature  were  chosen  until  the  value  of  temperature  for 
combustion  (Tc)  and  temperature  for  extinction  (Tg)  were 
approximately  equal  (Tc  «  Te)  or  (TQ  =  Te  +  10)  (F)  . 

C.  RESULTS 

For  each  case,  the  results  were  obtained  in  both  numer- 
ical and  graphical  form.  The  behavior  of  carbon  temperature, 
oxygen  concentration  and  reaction  rate  for  combustion  and 
extinction  are  shown  as  surfaces  versus  time  and  position. 

The  numerical  results  from  the  program  show  the  varia- 
tion of  several  parameter  values  during  the  transient  anal- 
ysis. The  air  temperature  inside  the  porous  medium,  the 
griphite  temperature  and  the  oxygen  concentration  are 
presented  at  each  position  (X/L)  at  intervals  of  time.  Other 
system  properties  are  also  given  .  These  properties  charac- 
terize the  process  of  combustion  and  extinction  for 
different  thicknesses.  The  maximum  and  minimum  do  not  occur 
at  the  same  point  or  at  the  same  time. 
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The  graphical  results  are  given  in  the  form  of  surfaces. 
These  surfaces  show  the  behavior  of  carbon  temperature, 
oxygen  concentration  and  reaction  rate  versus  time  and  posi- 
tion (X/L)  for  the  combustion  and  extinction  processes. 

1.  CASE  III-l   Thickness  =0.25  inches 

The  values  of  the  initial  uniform  carbon  temperature 
bounding  the  extinction  and  combustion  processes  are: 

for  combustion  Tc  =  1100  F 
for  extinction  Te  =  1090  F 

The  development  with  time  of  the  behavior  of  the 
graphite  temperature,  oxygen  concentration  and  reaction  rate 
during  the  combustion  and  extinction  processes  are  shown  in 
Figures  3.1  to  3.6  The  transient  parametric  analysis  provide 
some  insight  into  the  characterers  of  this  particular  case. 
The  maximum  and  minimum  values  of  these  parameters  are  given 
in  Table  V  . 

2.  CASE  III-2   Thickness  =0.50  inches 

The  combustion  and  extinction  process,  in  this  case, 
start  with  the  following  uniform  carbon  temperatures: 

for  combustion  Tc  =  960  F 
for  extinction  Te  =  950  F 


For  this  case,  the  graphical  results  for  graphite 
temperature,  oxygen  concentration  and  reaction  rate  are 
shown  in  Figures  3.7  to  3.12  .  The  characteristics  of  the 
changing  parameters  in  this  case,  are  shown  by  the  maximum 
and  minimum  values  in  Table  VI  . 
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TEMPERATURE   SURFRCE   FROM   GRRF3E 


Figure    3.1        Temperature   vs    X/L    and    time 

for    thickness    =    0.25    in 

Initial    carbon    temperature    =    1100    F. 

TEMPERATURE   SURFACE   FROM   GRRF3E 


Figure    3.2         Temperature   vs    X/L   and    time 

for    thickness    =    0.25    in 

Initial    carbon   temperature    =    1090    F. 
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OXYGEN  CONC.  SURFRCE  FROM  GRRF3E 
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o 
5  .»- 

1    "' 
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• 

Figure  3.3    Oxygen  concentration  vs  X/L  and  time 

for  thickness  =  0.25  in 

Initial  carbon  temperature  =  1100  F. 

OXYGEN  CONC.  SURFRCE  FROM  GRRF3E 


Figure  3.4    Oxygen  concentration  vs  X/L  arid  time 

for  thickness  =  0.25  in 

Initial  carbon  temperature  =  1090  F. 
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REACTION  RATE  SURFACE  FROM  GRAF3E 


Figure"  3.5    Reaction  rate  vs  X/L  and  time 

for  thickness  =  0.25  in 

Initial  carbon  temperature  =  1100  F. 

REACTION  RATE  SURFACE  FROM  GRAF3E. 


Figure  3.6    Reaction  rate  vs  X/L  and  time 

for  thickness  =  0.25  in 

Initial  carbon  temperature  =  1090  F. 


54 


CO 

UJ 

X 

u 

2 

M 

C 

3 

r^ 

m 

E 

i— i 

CM 

c 

•H 

vO 

• 

o 

c 

• 

O 

•H 

•H 

r—l 

U 

E 

ii 

3 

to 

J 

C 

CO 

E 

3 

CM 

w 

o 

E 

<r 

2 

u 

•H 

oo 

*; 

X 

• 

u 

05 

f— i 

LO 

I— 1 

vO 

CN 

CM 

vO 

r^ 

vO 

eo 

CM 

O 

• 

o 

rH 

sD 

• 

as 

CM 

-J" 

1 

i— 1 

o 

CM 

CM 


CM 

co 

CM 

<r 

O 

O 

-J" 

m 

CM 

• 

r>* 

-X) 

CM 

i 

vO 

O 

nT 

Pi 
o 
u* 

>     CO 

CO 

W        LU 

J  u 
a  o 
<    pi 


c 
o 

•H 
4-1 

a 

c 


x 

03 


CO 


vO 


CM 


r^ 

LP) 

n 

vO 

I-** 

co 

CO 

O 

\£> 

O 

ro 

.—1 

• 

o 

CO 

\D 

T— 1 

CM 

i 

CM 

O 

CM 
CM 

CO 

r^. 

m 

lO 

m 

O 

r— 1 

m 

ON 

i— » 

. 

i—( 

o> 

un 

• 

*>o 

CM 

i 

CO 

O 

n 

^  XI 


to 

4J 
4-1 

OJ 

E 
cO 

u 

cd 


c 


cr 

1—S 

^— - 

rf 

■H 

w 

cr 

-^. 

U 

• 

w 

4-1 

4-1 

•H 

*J 

• 

r* 

4-1 

U-( 

H-l 

u 

01 

*w 

m 

y-i 

» — ' 

4-t 

•H 

01 

OJ 

^^ 

-a 

>*. 

X 

o 

>N 

X 

a 

4-1 

g 

<j 

4-1 

1— 1 

u 

•H 

3 

•H 

v~-' 

00 

o 

o 

£ 

o 

•H 

<D 

a 

I— 1 

in 

■H 

XI 

^ 

>-i 

at 

"3 

4J 

C3 

3 

3 

> 

i— < 

u 

CJ 

w 

w 

O 

V 

E 

in 

M 

0) 

C 

> 

*M 

<U 

01 

>-. 

?o 

c 

(U 

5-i 

u 

o 

<3J 

o 

a 

Cu 

a 

a. 

Pi 

u 

55 


TEKPERRTURE  5URPRCE  PROM  GRRFZE 


Figure  3.7    Temperature  vs  X/L  and  time 

for  thickness  =  0.50  in 

Initial  carbon  temperature  =  960  F. 

7lm=ERR7URE  SURFRCE  FRC^  GRRF3E 


Figure  3.8    Temperature  vs  X/L  and  time 

for  thickness  =  0.50  in 

Initial  carbon  temperature  =  950  F. 
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OXYGEN  CONC.  SURFACE  FROM  GRRF3E 


Figure  3 . 9 


Oxygen  concentration  vs  X/L  and  time 
for  thickness  =  0.50  in 
Initial  carbon  temperature  =  960  F. 


OXYGEN  CONC.  SURr=GE  rRCl  GRRF3E 


K  isrti 


Figure  3.10  ,  Oxygen  concentration  vs  X/L  and  time 
for  thickness  =  0.50  in  __ 

Initial  carbon  temperature  =  950  F. 
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REACTION  RRTE  SURFf,C£.  ^  ^^ 


^gure  3.11    React.cn  rat, 


for  thickness  =  0.50  in 


'L  and  time 
initial  carbon  temperature  =  960  F. 
RATION  RATE  SURrfiCE  FROM  GRflr3E 


Figure  3.  12    Rearfl- 

Reactxon  rate  vs  X/L  and  time 

tor  thickness  =  0.50  in 

Initial  carbon  temperature  =  950  F. 
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CASE  III-3   Thickness 


0.75  inches 


During  this  case  the  minimum  uniform  temperature 
that  starts  the  combustion  process  and  the  maximum  uniform 
temperature  for  extinction  are: 

for  combustion  Tc  =  900  F 
for  extinction  Te  =  890  F 

During  the  transient  analysis,  the  system  parameters 
varied.  For  this  case  ,  the  maximum  and  minimum  values  of 
some  of  these  parameters  provide  some  insight  into  the  char- 
acter of  the  case.  These  parameters  and  their  extreme  values 
are  given  in  Table  VII  . 

Figures  3.13  to  3.18  show  the  graphical  results  for 
carbon  temperature,  oxygen  concentration  and  reaction  rate 
for  this  case. 

TEMPERRTURE  SURFACE  rROn  GRRF3E 


Figure  3.13    Temperature  vs  X/L  and  time 

for  thickness  =  0.75  in 

Initial  carbon  temperature  =  900  F. 
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TEMPERATURE  SURFACE  PROM  GRAF3E 


Figure  3.14    Temperature  vs  X/L  ana  time 

for  thickness  =  0.75  in 

Initial  carbon  temperature  =  890  F. 

OXYGEN  CONC.  SURFACE  FROM  GRAF3E 


Figure  3.15    Oxygen  concentration  vs  X/L  and  time 

for  thickness  =  0.75  in 

Initial  carbon  temperature  =  900  F. 
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OXYGEN  CONC.  SURFACE  FROM  GRAF3E 


Figure  3 . 16 


Oxygen  concentration  vs  X/L  and  time 
for  thickness  =  0.75  in 
Initial  carbon  temperature  =  890  F. 


REACTION  RATE  SURFACE  FROM  GRAF3E 


Figure  3.17    Reaction  rate  vs  X/L  and  time 

for  thickness  =  0.75  in 

Initial  carbon  temperature  =  900  F. 
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REACTION  RATE  SURFACE  FRCM  GRAF3C 


Figure  3.18    Reaction  rate  vs  X/L  and  time 

for  thickness  =  0.75  in 

Initial  carbon  temperature  =  890  F. 

4.   CASE  III-4   Thickness  =  1.00  inches 

The  program  was  run  with  several  values  of  initial 
uniform  graphite  temperature  .  The  temperatures  bounding 
combustion  and  extinction  are: 


for  combustion  Tc  =  850  F 
for  extinction  Te  =  840  F 

The  variable  parameters  of  this  case  attain  maximum 
and  minimum  values  as  given  in  Table  VIII  . 

In  this  case,  the  graphical  results  for  carbon 
temperature,  oxygen  concentration  and  reaction  rate  are 
shown  in  Figures  3.19  to  3.24  . 
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TEMPERATURE  SURFACE  FRCr  GRR~3E 


l:nc  i 


SCC; 


Figure  3.19    Temperature  vs  X/L  and  time 

for  thickness  =  1.00  in 

Initial  carbon  temperature  =  850  F. 


TEMPERATURE  SURFACE  FRCiM  GRA~3E 


Figure  3.20    Temperature  vs  X/L  and  time 

for  thickness  =  1.00  in 

Initial  carbon  temperature  =  840  F. 
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OXYGEN   CCNC.    SURrflCE   PROM   GRPE3E 


Figure    3.21         Oxygen    concentration  vs    X/L    and    time 

for    thickness    =    1.00    in 

Initial    carbon   temperature    =    850    F. 

:xygen  conc.  surface  erom  gpreze 


Figure  3.22    Oxygen  concentration  vs  X/L  and  time 

for  thickness  =  1.00  in 

Initial  carbon  temperature  =  840  F. 
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REACTION   RATE   SURFACE   FROH   GRRF3E 
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Figure    3.23         Reaction    rate   vs    X/L    and    time 

for    thickness    =    1.00    in 

Initial    carbon    temperature  '=    850    F. 

REACTION   RATE   SURrPCE   FROM   GRA~3E 
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Figure  3.24    Reaction  rate  vs  X/L  and  time 

for  thickness  =  1.00  in 

Initial  carbon  temperature  =  840  F. 
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5.   CASE  III-5   Thickness 


2.00  inches 


The  combustion  and  extinction  process,  in  this  case, 
are  characterized  by  the  following  uniform  initial  graphite 
temperatures : 

for  combustion  Tc  =  760  F 
for  extinction  Te  =  750  F 

The  graphite  temperature  and  other  parameters  varied 
during  these  processes  .  The  variation  of  these  parameters 
characterize  the  system.  The  maximum  and  minimum  values  of 
these  parameters  are  in  Table  IX  . 

The  graphical  results  for  carbon  temperature,  oxygen 
concentration  and  reaction  rate  are  shown  in  Figures  3.25  to 
3.30  . 

TEMPERATURE  SURFRCE  PROM  GRRP3E 


Figure  3.25    Temperature  vs  X/L  and  time 

for  thickness  =  2.00  in 

Initial  carbon  temperature  =  760  F. 
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TEMPERATURE  SURFACE  FROM  GRAF3E 
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Figure  3.26    Temperature  vs  X/L  and  time 


for  thickness  =  2.00  in 
Initial  carbon  temperature  =  750  F 

OXYGEN  CONC.  SURFPCE  FROM  GRAF3L 


Figure  3.27    Oxygen  concentration  vs  X/L  and  time 

for  thickness  =  2.00  in 

Initial  carbon  temperature  =  760  F. 
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OXYGEN   CONC.    SURFACE   FROM   GRRF3E 


Figure    3.28         Oxygen    concentration   vs    X/L    and    time 

for    thickness    =    2.00    in 

Initial    carbon   temperature    =    750    F. 

RERCTION  RRTE  SURFACE  FROM  GRfiiq££ 


Figure    3.29         Reaction   rate    vs    X/L    and    time 

for    thickness    =    2.00    in 

Initial    carbon   temperature    =    760    F. 
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REACTION  RATE  SURFACE  PROM  GRAP3E 
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Figure  3.30    Reaction  rate  vs  X/L  and  time 

for  thickness  =  2.00  in 

Initial  carbon  temperature  =  750  F. 

6 .   CASE  III-6   Thickness  =    4.00  inches 

The  minimum  value  of  the  initial  graphite  tempera- 
ture to  start  combustion  of  the  porous  medium  and  the 
maximum  value  of  the  initial  graphite  temperature  that 
results  in  extinction  for  this  case  are: 


for  combustion  Tc  =  690  F 
for  extinction  Te  -  680  F 

Other  variable  parameters  of  this  case  are  shown  in 
Table  X  .  The  values  of  these  parameters  are  characteristic 
of  the  combustion  and  extinction  processes  for  this  value  of 
thickness . 

In  this  case,  the  graphical  results  show  the 
behavior  of  carbon  temperature  ,  oxygen  concentration  and 
reaction  rate  in  Figures  3.31  to  3.36 
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"DMPE&FTu^E    SURfACE    FROn   GRRr3E 


J      N 


<tr<* 


Figure  3.31    Temperature  vs  X/L  and  time 

for  thickness  =  4.00  in 

Initial  carbon  temperature  =  690  F. 


'"3rs3- 


rL   5Crit  RCE   FRC~   GR»~-*r 


*  -    'A  '  V  y-.AXvt./.xVx 


Figure  3.32    Temperature  vs  X/L  and  time 

for  thickness  =  4.00  in 

Initial  carbon  temperature  =  680  F. 
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3£n  ::'.c.  suprpcc  from  s^^riE 


Fieure  3.33    Oxygen  concentration  vs  X/L  and  time 


for  thickness  =  4.00  in 
Initial  carbon  temperature  =  690  F 


D'/-3e:n  :cr,c.  susrrCE  from  3RRr3r_ 


Figure  3.34    Oxygen  concentration  vs  X/L  and  time 


for  thickness  =  4.00  in 
Initial  carbon  temperature  =  680  F 
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:erctig\  PR?:  3u=r-::  pro.-  &aafze 


Figure  3.35    Reaction  rate  vs  X/L  and  time 
for  thickness  =  4.00  in- 
Initial,  carbon  temperature  .-  690  F. 


REACTION  RATE  SURFACE  FROM  C-RRF3E 


Figure  3.36    Reaction  rate  vs  X/L  and  time 

for  thickness  =  4.00  in 

Initial  carbon  temperature  =  680  F. 
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CASE  III- 7   Thickness 


6 . 00  inches 


Below  is  given  the  initial  uniform  carbon  tempera 
ture  for  combustion  and  extinction: 


for  combustion  Tc  =  650  seconds 
for  extinction  Te  =  640  seconds 

The  variable  parameters  for  this  case  are  given  in 
Table   XI 

The  Figures  3.37  to  3.42  show  the  development  of 
carbon  temperature,  oxygen  concentration  and  reaction  rate 
at  each  position  during  the  combustion  and  extinction 
processes . 


Figure  3.37    Temperature  vs  X/L  and  time 

for  thickness  =  6.00  in 

Initial  carbon  temperature  =  650  F. 
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TEMPERATURE  SURFRCE  FROM  GRRF3E 


Figure  3.38    Temperature  vs  X/L  and  time 


for  thickness  =  6.00  in 
Initial  carbon  temperature  =  640  F 


Figure  3.39    Oxygen  concentration  vs  X/L  and  time 

for  thickness  =  6.00  in 

Initial  carbon  temperature  =  650  F. 
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OXYGEN  CONC.  SURFACE  FROM  GRRF3E 


Figure  3.40    Oxygen  concentration  vs  X/L  and  time 

for  thickness  =  6.00  in 

Initial  carbon  temperature  =  640  F. 


Figure  3.41    Reaction  rate  vs  X/L  and  time 

for  thickness  =  6.00  in 

Initial  carbon  temperature  =  650  F. 
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REACTION   RRTE   SURFRCE  FRCM   G^flFZE 


Figure    3.42         Reaction    rate   vs    X/L   and    time 

for    thickness    =    6.00    in 

Initial    carbon    temperature    =    640    F. 

D .       SUMMARY 

Here,  some  observations  are  made  about  the  effects  of 
thickness  on  the  combustion  process. 
1 .   Power  Relation 

For  each  case  a  pair  of  values  (thickness,  combus- 
tion temperature)  was  obtained.  Four  plots  of  this  data  are 
shown  in  Figures  3.43  to  3.46  . 

Figures  3.43  and  3.44  show  thickness  versus  combus- 
tion temperature  ,  and  Figures  3.45  and  3.46  show  thickness 
versus  the  time  it  takes  for  the  oxygen  concentration  at  X/L 
=  1.  to  reach. 

It  is  observed  that  the  relation  between  thickness 
and  combustion  temperature  ,  figure  3.43  has  a  form  which 
approximates  an  exponentially  decreasing  function. 

The  log-log  plot  (  Figure  3.44  )  is  close  to  a 
linear  form,  which  yields  the  approximate  power  relation. 

The  temperature  is  in  degrees  fahrenheit  and  thick- 
ness is  in  inches  . 

81 


c/j 

X 

u 

2 


O 
O 


c 
o 

■H 

j-j 

CO 

3 
XI 
c 

o 

u 


3 
E 
■H 
C 


c 

3 

s 

•H 

CO 

e 


^o 


CM 


00 


vi5 


• 

m 

vD 

ON 

o 

^O 

<f 

o 

O 

r-l 

CN 

• 

CM 

I 

ON 

o 

vO 


<r 

O 

r^ 

• 

CM 

en 

m 

CM 

O 

. 

m 

ON 

• 

o 

i 

CM 

o 

r— 1 

3 

3 

3     -h 
O       3 

•H       -H 

xj  £ 
u 
3 

•H  3 

xj  3 

X!  £ 

CD       -H 
X 

£ 


CO 


vO 


r^ 

O 

vO 

vO 

00 

o 

r— 1 

m 

CM 

1 

i— i 

CN 


CO 


CO 


CN 


CM 


CO 


en 

o 


r>» 


w 

s-i 

XJ 

<u 

e 

c5 


— 1 

■>■»» 

XJ 

' — ' 

X 

- — v 

3 

' >> 

i— i 

u 

<D 

CT 

<*"N 

^^ 

x 

■H 

w 

CT 

■-«, 

U 

m 

XJ 

XJ 

•H 

XJ 

3 

<4-i 

CH 

4-1 

XJ 

OJ 

■> — ' 

u 

(XJ 

^— ' 

<XJ 

■H 

01 

OJ 

"•■», 

-a 

>N 

X 

o 

>> 

X 

a 

XJ 

£ 

u 

XJ 

T— 1 

u 

•H 

3 

•H 

» — ' 

Qfl 

O 

3 

r* 

r— 1 

O 

o 

•H 

OJ 

(U 

i— ( 

CO 

•H 

XI 

5-i 

}-i 

<v 

T3 

XJ 

rt 

3 

3 

> 

r-H 

<J 

oj 

to 

CO 

O 

CD 

c 

w 

CO 

0J 

3 

> 

U 

CD 

0J 

i-l 

^ 

3 

CD 

U 

i-l 

o 

cd 

o 

a 

a 

a 

a 

Pi 

CJ 

82 


T  =  850  L-0.16554  (31) 

A  comparison  of  the  computer  program's  results  and 
the  formula's  results  are  given  in  Table  III.  The  maximum 
difference  between  the  analytical  results  and  the  power 
equation  results  is  less  than  3  percent  over  the  range  of 
thickness  between  0.25  and  6.00  inches  . 

This  shows  that  the  formula  is  a  good  approximation 
for  these  cases. 

2 .   Combustion  Start  Time 

During  the  initial  phase  of  the  combustion  process, 
the  oxygen  concentration  at  all  points  of  the  porous  medium 
continuously  decreases  to  zero,  for  all  cases,  the  point  of 
the  porous  medium  that  first  reached  zero  was  X/L  =  1.  The 
time  it  taken  for  that  the  oxygen  concentration  at  X/L  =  1. 
to  reach  zero  is  important  because  it  is  related  to  the  time 
that  combustion  starts. 

Figure  3.45  ,  shows  that  this  time  does  not  increase 
linearly  with  the  thickness.  The  log- log  plot  of  these 
values,  Figure  3.46  ,  shows  that  this  log-log  plot  is 
approximately  a  straight  line. 

The  equation  for  this  line  is: 

T*  =  226.8  L1.683  (3.2) 

This  equation  can  be  used  to  give  an  estimate  for 
the  time  it  takes  for  combustion  to  start. 
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THICKNESS  VERSUS  TEMPERATURE 


8       3       4       5 

THICKNESS  (INCHES) 


Figure  3 . 43 


Thickness  versus  combustion  temperature 
Retangular  plot. 
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THICKNESS  X  TEMPERATURE 


THICKNESS  (INCHES) 


Figure  3.44    Thickness  versus  combustion  temperature 

Log- log  plot . 
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THICKNESS  VERSUS  T* 


t       9       4        ft 

THICKNESS  (INCHES) 


Figure  3 . 45 


Thickness  vs  time  for  oxygen  concentration 
at  X/L  =  1.  to  reach  zero 
Retangular  plot. 
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THICKNESS  VERSUS  T* 


THICKNESS  (INCHES) 


Figure  3.46    Thickness  vs  time  for  oxygen  concentration 

reach  zero  log-log  plot. 


87 


oo 


vo     m     oo     rn     vo 


CM 

CO 

m 

CM 

<t 

m 

CM 

CNJ 

vO 

i— i 

CM 

r>^ 

cm 
CM 

CO 

5-i 

w 

O 

C^ 

O 

m 

r-l 

CO 

CO 

5-i 

CO 

r>» 

<! 

5-4 

• 

O 

O 

O 

CM 

CM 

u 

[n 

0) 

CM 

O 

CO 

* — ' 

J-S 

co 

u 

01 

z 

u 

^ 

3 

1-1 

u 

U 

M 

tH 

rt 

X, 

— 

5-1 

5-i 

H 

<U 

<U 

£J 

o 

O 

o 

O 

O 

w 

a 

5 

o 

O 

O 

ON 

<r 

m 

00 

<r 

J 

J 

g 

o 

•H 

ON 

m 

CO 

CO 

I**. 

\o 

o 

cq 

hJ 

3 

G. 

4J 

o 

o^ 

< 

< 

XJ 

cd 

t—l 

H 

o 

CO 

o 

■H 

J-J 
W 

i— 1 

S-i 

g 

H 

3 

rt 

o 

o 

o 

o 

o 

o 

o 

J 

J3 

5-i 

o 

vO 

o 

m 

o 

o^ 

in 

3 

S 

aO 

1—< 

o> 

o> 

CO 

r> 

\£> 

vO 

CO 

O 

o 

1— 1 

u 

u 

5-i 

c* 

C 

o 

•H 

CU 
5-i 

4-1 

a3 

a 

U 

5-( 

o 

O 

o 

o 

o 

o 

O 

G 

OJ 

' — v 

o> 

U~> 

a> 

^~ 

m 

oo 

<r 

■H 

a 

ft, 

o 

ON 

00 

00 

r>* 

o 

o 

JJ 

s 

» — ' 

T—l 

X 

01 

01 

4J 

c 

m 

o 

\T\ 

o 

O 

o 

o 

•H 

CM 

m 

r^ 

o 

O 

o 

o 

J 

o 

o 

o 

I—I 

CM 

<r 

\D 

w 

1— 1 

CM 

m 

<r 

m 

\D 

r^. 

OJ 

1 

l 

i 

i 

i 

1 

i 

CO 

M 

H 

i-i 

H 

i— i 

1— 1 

H 

Kl 

H 

H 

H 

H 

H 

H 

H 

U 

M 

M 

M 

H 

H 

H 

l-l 

88 


IV.  PERMEABILITY 

A.   INTRODUCTION 

The  permeability  or  hydraulic  conductivity  of  a  porous 
medium  is  the  capacity  of  air  flow  through  the  porous 
medium.  The  air  provides  the  oxygen  for  combustion  .  On  the 
other  hand  ,  air  flow  provides  for  convective  heat  transfer 
within  the  medium  .  Thus  air  flow  provides  for  heat  genera- 
tion from  combustion  ,  as  well  as  heat  transfer  by  convec- 
tion .  The  interaction  of  these  two  mechanisms  determines 
whether  combustion  or  extinction  will  occur  .  If  the  heat 
generation  of  combustion  is  greater  than  the  heat  transfer 
of  convection  ,  then  combustion  will  occur  .  If  the  heat 
transfer  mechanism  dominates  ,  then  extinction  will  occur  . 
The  equation  used  for  permeability  in  the  program  is  : 

m  -  0.2  p3  /  z2  (4.1) 

Where  m  is  the  permeability  ,  p  is  the  porosity  and  z  is  the 
specific  internal  area  .  The  expressions  for  porosity  and 
specific  internal  area  are  : 

p  =  1  -  v    I    6(d/D)3  (4.2) 


z  =  1/2  77-  (d/D2)  (4.3) 


z  =  1/2  v    (d2/D3)  (4.4) 

The  equation  4.3  is  for  spherical  particles  and  equation  4.4 
is  for  cylindrical  particles  .   Equation  4.1  is  the  Kozeny 
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equation  for  permeability  of  porous  medium  .  This  equation 
fails  for  highly  porous  fibrous  media  . 


Figure  4.1    Geometry  of  a  typical  cell. 

This  equation  fails  for  highly  porous  fibrous  medias  . 

The  tortuosity  ,t  ,  is  the  ratio  of  the  length  of  the 
flow  path  for  a  fluid  particle  to  the  straight  line 
distance.  The  tortuosity  depends  on  the  ratio  d/D  .In  this 
program  tortuosity  was  assumed  to  have  a  constant  value  of 
1.4   . 

Normally  the  interior  geometry  of  a  porous  medium  can  be 
inregular  and  complex  .  In  such  cases  the  internal  structure 
is  difficult  to  discribe  analytically  .  The  program  used  in 
this  investigation  models  the  porous  medium  as  consisting  of 
cylindrical  fibers  or  spherical  particles  in  a  regular  peri- 
odxc  geometry  as  shown  in  Figures  4.1  . 

During  the  combustion  process  carbon  is  being  continu- 
ously consumed  and  as  a  result,  the  internal  geometry  is 
changing.  The  medium  becomes  more  porous  as  the  carbon 
decreases.  Therefore,  all  properties  which  depend  upon  fiber 
or  particle  diameter  are  functions  of  time  and  position.  The 
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model  used  in  the  program  assumed  that  the  carbon  matrix 
remains  rigid  as  the  particles  diameter  decreases  ,  and  thus 
,  the  porosity  increases  with  combustion  . 

B .   PROCEDURE 

To  study  the  effects  of  permeability  on  combustion  ,  all 
parameters  of  the  program  were  maintained  constant  ,  except 
the  permeability  .  In  order  to  change  the  permeability  and 
keep  the  porosity  constant  ,  it  is  necessary  to  keep  the 
(d/D)  ratio  constant  as  seen  in  the  equations  4.5  and  4.6  . 
This  is  accomplished  by  setting  d/D  equal  too:  .  Then  the 
numerator  of  equation  4.1  remains  constant  ,  however  the 
denominator  changes  because   d  =a   D  and  Z  gives  : 

Z  =  1/2  7r(a  d/D2)  =  constant  /  D  (4.5) 


Z  =  1/2  <n{a2    d2/D3)=  constant'/  D  (4.6) 

Thus  ,  the  dimensions  of  d  and  D  were  changed  propor- 
tionally so  that  the  initial  ratio  (d/D)  remained  constant 
in  all  cases  .  Thus  the  porosity  remained  constant  while  the 
permeability  was  different  for  each  case  . 

The  cases  studied  in  this  section  are  : 

CASE  IV-1  PERMEABILITY  =  0.00015  ft2  (d=D=0.0025  in) 

CASE  IV-2  PERMEABILITY  =  0.00058  ft2  (d=D=0.0050  in) 

CASE  IV-3  PERMEABILITY  =  0.00233  ft2  (d=D=0.0100  in) 

CASE  IV-4  PERMEABILITY  =  0.00933  ft2  (d=D=0.0200  in) 

For  each  values  of  permeability  the  program  was  run  with 
different  uniform  carbon  temperatures  for  all  points  (X/L)  . 
This  procedure  was  repeated  until  we  obtaind  the  minimum 
temperature  that  results  in  combustion  .   This  temperature 
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is   called  the   combustion  temperature  (Tc)   .   Extinction 
occurs  at  temperatures  below  T   . 

Tc  =  Te  +  S  (4.7) 

In  this  section  this  difference  (<5  )  was  chosen  equal  at 
10  F  then  : 

Tc  =  Te  +  10   (F)  (4.8) 


C.   RESULTS 

For  each  case  ,  the  program  provided  the  values  of  the 
carbon  and  air  temperatures  inside  the  porous  medium  ,  as 
well  as  other  parameters  which  varied  during  the  combustion 
process  . 

1.   CASE  IV- 1  Permeability  =  0.00015  ft2  (d=D=0.0025  in) 

For  this  value  of  permeability  ,  the  maximum  carbon 

temperature  that  does  not  result  in  combustion  (extinction 

temperature)  and  the  minimum  temperature  that  results  in 

combustion  (combustion  temperature)  are  : 


for  combustion     T   =  760   F 
for  extinction     Te  =  750   F 


The  graphical  results  are  shown  in  Figures  4.2  to 
4.7  .  These  graphical  surfaces  show  the  carbon  temperature  , 
oxygen  concentration  and  reaction  rate  at  each  position  and 
time  .  The  form  of  these  surfaces  shows  the  different 
behavior  of  these  three  parameters  during  the  extinction  and 
combustion  processes. 
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TEMPERATURE  SURFACE   FROM   GRAF3E 


Figure    4.2         Temperature   vs    X/L    and    time 

for   permeability    =    0.00015    ft2 

initial    carbon    temperatiure    =    760    F. 


TEMPERATURE   SURFACE   FROM   GRPF3E 


Figure    4.3         Temperature   vs    X/L    and    time 

for   permeability    =    0.00015    ft2 

initial    carbon    temperatiure    =    750    F. 
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OXYGEN   CONC.    SURFACE   FRCH   GRAF3E 


Figure    4.4         Oxygen    concentration   vs    X/L   and    time 

for   permeability    =    0.00015    ft2 

initial    carbon   temperatiure    =    760    F. 


OXYGEN   CONC.    SURFACE   FROM    GRPF3E 


Figure    4.5         Oxygen    comcentrat ion   vs    X/L   and    time 

for   permeability    =    0.00015    ft2 

initial    carbon    temperatiure    =    750    F. 
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REACTION   RATE   SURFACE   FROM   GRPF3E 
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Figure    4.6         Reaction   rate   vs    X/L    and    time 
for   permeability    =    0.00015    ft 
initial    carbon   temperatiure    =    760    F. 


REACTION   RATE   SURFACE   FROM   GRPF3E 


Figure    4.7         Reaction    rate   vs    X/L    and    time 
for   permeability    =    0.00015    ft 
initial    carbon   temperatiure    =    750    F. 
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2.   CASE  IV-2  Permeability  =  0.00058  ftz  (d=D=0.0050  in) 

n 
Here  the  value  of  permeability  is  0.00058  ft   which 

results  from  a  pore  diameter  of  0.0050  in  .  For  this  case 

the  combustion  and  extinction  temperatures  are  : 


for  combustion 
for  extinction 


Tc  =  850  0  F 
Te  =  840  0  F 


From  the  numerical  output  of  the  program  ,  three 
graphical  surface  s  were  created  ;  the  carbon  temperature  , 
oxygen  concentration  and  reaction  rate  as  functions  of  time 
and  position  .  These  surfaces  provide  information  of  the 
development  of  these  three  dependent  variables  during  the 
combustion  and  extinction  process  .  These  surfaces  are  shown 
on  Figures  4.8  to  4.13  . 


TEMPERATURE  SURERCE  ERCM  GRRP3E 


Figure  4.8    Temperature  vs  X/L  and  time 

for  permeability  =  0.00058  ft2 

initial  carbon  temperatiure  =  850  F. 
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TEMPERATURE   SURFACE   FROM   GRP^3E 
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Figure    4.9         Temperature   vs    X/L    and    time 

for   permeability    =    0.00058    ft2 

initial    carbon    temperatiure    =    840    F. 


OxrGEN   CONC.    SURFACE   EPOn   GRAF3E 


Figure    4.10         Oxygen    concentration  vs    X/L    and    time 

for   permeability    =    0.00058    ft2 

initial    carbon   temperatiure    =    850    F. 
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OXYGEN   CONG.    SURFRCE    FROM   GRRF3E 


Figure    4.11         Oxygen    comcentration   vs    X/L    and    time 

for   permeability    =    0.00058    ft2 

initial    carbon    temperatiure    =    840    F. 


PERCTI0N   RRTE    SURFRCE   FROM    GRRr3E 


Figure    4.12         Reaction   rate   vs    X/L   and    time 

for   permeability    =    0.00058    ft2 

initial    carbon    temperatiure    =    850    F. 
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REACTION  RATE  SURFACE  PROM  GRAF3E 
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Figure  4.13    Reaction  rate  vs  X/L  and  time 

for  permeability  =  0.00058  ft2 

initial  carbon  temperatiure  =  840  F. 

3.   CASE  IV-3  Permeability  =  0.00233  ft2  (d'=D=0.0100  in) 
For  this  case  the  initial  permeability  is  0.00233 
ft   .  With  this  value  of  permeability  the  combustion  and 
extinction  temperatures  are  : 


for  combustion 
for  extinction 


Tc  =  980  0.F 
Te  =  970  0.F 


The  variation  of  the  dependent  variables  during  the 
combustion  and  extinction  cases  is  shown  by  a  numerical 
output  of  the  program  .  These  results  are  also  shown  in 
three  graphical  surfaces  ;  carbon  temperature  ,  oxygen 
concentration  and  reaction  rate  on  each  position  (X/L)  and 
time  .  These  graphical  results  are  shown  in  Figures  4.14  to 
4.19 
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TEMPERATURE  5URFRCE   ^RCr   3RRF3E 


Figure    4.14         Temperature   vs    X/L    and    time 

for   permeability    =    0.00233    ft2 

initial    carbon    temperatiure    =    980    F. 

TEMPERATURE   5URFPCE   FRCT   GRPP3E 


Figure    4.15         Temperature    vs    X/L    and    time 

for   permeability    =    0.00233    ft2 

initial    carbon    temperatiure    =    970    F. 
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OXtGEN  CONC.  SURFRCE  rRC^  GRQF3E 
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Figure  4.16    Oxygen  concentration  vs  X/L  and  time 

for  permeability  =  0.00233  ft2 

initial  carbon  temperatiure  =  980  F. 


OXYGEN  CONC.  SURFACE  rRC^  G*CF3E 


tT*  '*c' 


Figure  4.17    Oxygen  comcentrat ion  vs  X/L  and  time 

for  permeability  =  0.00233  ft2 

initial  carbon  temperatiure  =  970  F. 
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^lRCTION  RATE  5UR~pCE  FRCM  G3Rr3E 


Fieure  4.18    Reaction  rate  vs  X/L  and  time 

2 
for  permeability  =  0.00233  ft 

initial  carbon  temperatiure  =  980  F. 


REACTION  RRTE  SURFRCE  ERCf  GRRr3i 


Figure  4.19    Reaction  rate  vs  X/L  and  time 
for  permeability  =  0.00233  ft 
initial  carbon  temperatiure  =  970  F. 
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4.   CASE  IV-4  Permeability  =  0.00933  ft2  (d=D=0.0200  in) 

The  combustion  and  extinction  temperatures  for  this 

o 
case  ,  with  initial  permeability  equal  to  0.00933  ft   ,  are: 


for  combustion    Tc  =  1180  0.F 
for  extinction    Te  =  1170  0.F 


for  this  case  ,  it  was  not  possible  to  obtain  the  graphical 
surfaces  due  to  limitations  of  the  graphical  program  . 

D .   SUMMARY 

Here  ,  some  observations  about  the  effects  of  perme- 
ability on  combustion  can  be  understood  by  looking  at  all 
the  results  at  one  time  .  Table  XIII  contains  the  results  of 
all  cases  . 

The  Table  lists  the  temperatures  for  extinction  and 
combustion  ,the  permeability  (m)  ,  the  dimension  of  the  unit 
cell  (d  =  D)  and  time  t*  .  The  t*  time  is  the  time  at  which 
the  oxygen  concentration  at  a  point  reaches  zero  first  . 
That  point  of  the  porous  medium  that  first  reaches  zero 
oxygen  concentration  in  all  these  cases  is  the  the  position 
point  X/L  =  1  . 

1 .   Power  Relation 

The  pair  of  values  (permeability  ,  combustion 
temperature)  for  each  case  was  plotted  on  a  cartesian  graph. 
An  observation  of  this  curve  shows  that  the  combustion 
temperature  increases  with  increasing  permeability.  Figure 
4.20  shows  that  the  combustion  temperature  increases  mono- 
tonically  and  approaches  an  assindot icptate  as  permeability 
increases  .  The  plot  of  these  same  points  on  a  log- log 
graph  ,  Figures  4.21  ,  shows  that  its  form  is  not  quite  a 
straight  line  . 

It  is  observed  that  the  relation  between  combustion 
temperature  (T  )  and  permeability  (m)  yields  the  approximate 
power  relation  (equation  4.9)  : 
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Tc  =  1880  m0.105  (4.9) 

A  comparation  of  the  values  of  the  combustion 
temperature  obtained  by  the  program  and  by  the  power  rela- 
tion is  presented  in  Table  XIV  .  This  Table  shows  that  the 
difference  between  these  results  fail  between  1.2  to  3.2 
percent  .  This  result  shows  that  the  above  power  relation  is 
a  good  approximation  to  estimate  the  temperature  which  will 
start  combustion  for  a  given  permeability  . 
2 .   Combustion  Speed 

In  this  section  ,  the  question  of  how  permeability 
affects  the  time  it  takes  for  combustion  to  begin  is 
addressed  .  We  meassure  the  rapidilly  of  the  combustion 
process  by  the  t*  time  previously  defined  .  Recall  that  t* 
is  the  time  at  which  the  oxygen  first  reaches  zero  at  some 
point  in  the  porous  medium  .  In  all  cases  the  position  X/L  = 
1  is  the  first  point  to  achieve  zero  oxygen  concentration  . 
The  plots  of  permeability  (m)  versus  t-  ,  Figures  4.22  and 
4.23  ,  show  that  t*  decays  exponentially  with  increasing 
permeability  and  approaches  an  assymptotic  value  of  t-  =  0 
as  m  goes  to  infinity  . 

A  very  interesting  observation  is  obtained  from  the 
results  of  Table  XIII  .  It  is  noted  that  the  m^Wmj  (for  i 
=  1,2,3)  remains  pratically  constant  (3.84  ,  4.01  ,  4.00) 
while  the  t'"-+i/t*-  ratio  takes  on  the  variable  values 
l/(i+l)  (1/1.91  ,  1/2.65  ,  1/3.56)  . 

An  attempt  was  made  to  obtain  an  additional  run  with 
the  permeability  equal  to  0.03733  ft  ,  however  the  combus- 
tion activity  was  too  rapid  to  achieve  meaningful  results  . 
Combustion  was  achieved  so  rapidly  that  the  gradient  were 
too  large  for  the  program  to  handle  . 

The  permeability  is  a  parameter  that  has  great 
influence  on  the  combustion  process  for  porous  medium  . 
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PERMEABILITY  VERSUS  TEMPERATURE 


0.000  0.00*  0.004  0.OM  0.000 

PERMEABILITY  (FT.SQ) 
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Figure  4.20 


Permeability  versus  combustion  temperature 
retangular  plot . 
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PERMEABILITY  VERSUS  TEMPERATURE 


PERMEABILITY  (FT.SQ) 


Figure  4.21    Permeability  versus  combustion  temperature 


log  log  plot. 
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PERMEABILITY  VERSUS  T* 


0.000 


0.001  0.004  0.006  0. 

PERMEABILITY  (FT.SQ) 


0.010 


Figure  4.22    Permeability  versus  t' 
Retangular  plot. 
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PERMEABILITY  VERSUS  T* 


PERMEABILITY  (FT.SQ) 


Figure    4.23         Permeability   versus    t* 
Log-log    plot. 
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V.  REACTION  ORDER 

A.   INTRODUCTION 

The  expression  for  reaction  rate  for  carbon  reacting  in 
air  used  in  the  program  was  the  expression  resulting  from 
the  experimental  data  of  Parker  and  Hotel  . 


gr-c 


rr  (5.1) 


Cnrvs 


Reaction  order  is  the  exponent  of  the  partial  pressure 
of  oxygen  in  the  expression  of  the  reaction  rate  .  The 
program  permits  the  user  to  select  any  value  for  this  param- 
eter .  In  all  of  the  previous  chapters  ,  a  value  of  n  equal 
to  0.5  was  used  .  In  this  chapter  we  investigate  the  effect 
of  n  on  the  combustion  process  . 

B .   PROCEDURE 

In  this  section  ,  all  parameters  are  fixed  to  their 
basic  values  and  only  the  reaction  order  is  changed  in  each 
case  . 

The  program  was  run  ,  for  each  value  of  reaction  order  , 
with  different  value  of  the  uniform  initial  carbon  tempera- 
ture until  we  obtained  the  minimum  temperature  which  results 
in  combustion  .  This  temperature  is  the  combustion  tempera- 
ture (T  )  .  The  maximum  temperature  that  results  in  extinc- 
tion is  the  extinction  temperature  (T  )  .  they  differ  by  an 
inf initesimal5  (Tc  =  Te  +5)  •  The  magnitude  of<5  was  chosen 
equal  to  10  F  for  this  study  . 

In  order  to  study  the  influence  of  the  reaction  order  on 
the  combustion  process  ,  four  different  values  of  this 
parameter  were  chosen  .  These  cases  are  : 
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CASE  V-l  Reaction  Order  =0.25 

CASE  V-2  Reaction  Order  =0.50 

CASE  V-3  Reaction  Order  =  0.75 

CASE  V-4  Reaction  Order  =  1.00 

The  fixed  parameters  used  in  all  cases  are  : 
Ambient   temperature  =80   F 
Ambient  pressure  =  2117  lb/ft. sq 
Tortuosty   =  1.400 

Filamento  diameter   =  0.0004167  ft 
Thickness  of  matrix  laminette  =  0.0004167  ft 
Thickness  of  porous  medium  =  0.02083  ft 
Gas  constant  for  air  =  53.34  lbf  ft/lbm  R 
Conductivity  of  filament  =  86  BTU/lbm  H.F. 
Specific  heat  of  filament  =  0.703  lbm/Cf 
Emissivity  of  filament  =  0.90 
Shape  factor  for  internal  HFxFER  =1.00 
Characteristc  length  of  medium  =  1  inchess 
Heat  of  reaction   =  14090  BTU/Cf 
Stochiometric  ratio  (fuel/air)  =  0.375 
Reaction  coefficient  =  2065000  lbm/Cf  H 
Activetion  energy  coefficient  =  28840  deg  R 
Pressure  differential  across  thickness  =  50  lb/ft. sq 
Initial  uniform  oxygen  concentration   =  0.172  lbm/ft 

C.   RESULTS 

For  each  cases  the  results  were  obtained  in  both  numer- 
ical and  graphical  form  .  These  results  show  the  variation 
of  some  parameters  during  the  combustion  and  extinction 
processes  .  The  results  obtained  for  each  value  of  reaction 
order  are  presented  below 

1.   CASE  V-l   Reaction  Order  =0.25 

With  this  value  of  reaction  order  ,  the  combustion 
and  extinction  processes  define  the  following  characteristic 
temperatures  : 
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for  combustion   Tc  =  980  o.F 
for  extinction   Tg  =  970  o.F 


The  analysis  of  the  numerical  output  of  the  program 
shows  that  some  parameters  varied  during  the  combustion  and 
extinction  processes  .  The  limits  of  these  parameters  show 
the  character  of  the  particular  kind  of  process  which 
occured  . 

The  graphical  results  show  three  surfaces  that 
represent  the  evolution  of  the  carbon  temperature  ,  oxygen 
concentret ion  and  reaction  rate  at  each  position  during  the 
transient  history  .  These  graphical  surfaces  are  shown  in 
Figures  5.1  to  5.6  . 


TEMPERATURE  5URrPCE  F^dV   GRqr3E 


Figure  5.1    Temperature  vs  X/L  and  time 

for  reaction  order  =  0.25 

Initial  carbon  temperature  =  980  F. 
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TEMPERATURE  SURPRCE  PROM  GRRP3E 


Figure    5.2         Temperature   vs    X/L    and    time 

for    reaction   order    =    0.25 

Initial    carbon    temperature    =    970    F. 

DxrGEN   CONC.    SURFACE   rRC"   C-RQF3E 


t:-c  i sec i 


Figure  5.3    Oxygen  concentration  vs  X/L  and  time 

for  reaction  order  =  0.25 

Initial  carbon  temperature  =  980  F. 
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OXYGEN  CONC.  SURFACE  PROM  GRRF3E 
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Figure  5.4    Oxygen  concentration  vs  X/L  and  time 

for  reaction  order  =  0.25 

Initial  carbon  temperature  =  970  F. 

PERCTION  RRTE  SURFPCE  FPCM  GPRF3E 


Figure  5.5    Reaction  rate  vs  X/L  and  time 

for  reaction  order  =  0.25 

Initial  carbon  temperature  =  980  F. 
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REACTION  RATE  SURFACE  FROM  GRAF3E 


Figuie  5.6    Reaction  rate  vs  X/L  and  time 

for  reaction  order  =  0.25 

Initial  carbon  temperature  =  970  F. 

2.   CASE  V-2   Reaction  Order  =  0.50 


For  reaction  order  equal  to  0.50  ,  the  initial 
uniform  carbon  temperatures  that  separate  the  combustion  and 
extinction  processes  are  : 


for  combustion 
for  extinction 


Tc  =  980  F 
T„  =  970  F 


The  program's  numerical  output  included  the  values 
of  some  parameters  that  varied  during  the  transient  process. 
The  variation  of  the  values  characterize  the  combustion  or 
extinction  process.  The  graphical  rasults  are  shown  in 
three  surfaces:  carbon  temperature,  oxygen  concentration  and 
reaction  rate  versus  position  (X/L)  and  time.  These  graph- 
ical surfaces  are  shown  in  figures  5.7  to  5.12  . 
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TEMPERATURE   SURFACE   FROM   GRAF3E 
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Figure    5.7         Temperature   vs    X/L    and    time 

for    reaction    order    =    0.50 

Initial    carbon    temperature    =    850    F. 

TEMPERATURE   SURFACE   FROM   GRAF3E 
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Figure    5.8         Temperature   vs    X/L   and    time 

for   reaction   order    =    0.50 

Initial    carbon   temperature    =    840    F. 
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OXYGEN  CONC.  SURRRCE  PROM  GRRE3E 


Figure  5,9    Oxygen  concentration  vs  X/L  and  time 

for  reaction  order  =  0.50 

Initial  carbon  temperature  =  850  F. 

OXYGEN  CONC.  SURFRCE  FROM  GRPF3E 


Figure  5.10    Oxygen  concentration  vs  X/L  and  time 

for  reaction  order  =  0.50 

Initial  carbon  temperature  =  840  F. 
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REACTION   RATE   SURFACE   FROM   GRRF3E 


Figure    5.11         Reaction   rate   vs    X/L   and    time 

for    reaction   order    =    0.50 

Initial    carbon    temperature    =    850    F. 
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Figure    5.12         Reaction   rate   vs    X/L    and    time 

for    reaction   order    =    0.50 

Initial    carbon    temperature    =    840    F. 
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3.   CASE  V-3   Reaction  Order  =0.75 

The  reaction  order  was  set  equal  to  0.75  for  this 
case.  The  combustion  and  extinction  temperatures  are: 


for  combustion   Tc  =  750  F 
for  extinction   Te  =  740  F 

A  visualization  of  the  development  of  the  carbon 
temperature,  oxygen  concentration  and  reaction  rate  at  each 
position  and  time  are  shown  in  the  graphical  surfaces  of 
Figures  5 . 13  to  5 . 18  . 


TEMPERATURE  SURFACE  FROM  GRAF3E 


Figure  5.13    Temperature  vs  X/L  and  time 

for  reaction  order  =  0.75 

Initial  carbon  temperature  =  750  F. 


120 


TEMPERATURE   SURFACE    f~Rcr   GRAF3E 
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Figure    5.14         Temperature   vs    X/L    and    time 

for   reaction   order    =    0.75 

Initial    carbon    temperature    =    740    F. 

OXYGEN   CONC.    SURFACE   FRQr   GRAF3E 


Figure    5.15         Oxygen    concentration  vs    X/L    and    time 

for    reaction    order    =    0.75 

Initial    carbon    temperature    =    750    F. 
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OXYGEN  CONC.  SURFACE  FROM  GRAF3E 


Figure  5.16    Oxygen  concentration  vs  X/L  and  time 


for  reaction  order  =  0.75 
Initial  carbon  temperature  -    740  F 

REACTION  RATE  SURFACE  FROM  GRAF3E 


Figure  5.17    Reaction  rate  vs  X/L  and  time 

for  reaction  order  =  0.75 

Initial  carbon  temperature  =  750  F. 
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REACTION  RATE  SURFACE  FROM  GRAF3E 


Figure  5.18    Reaction  rate  vs  X/L  and  time 

for  reaction  order  =  0.75 

Initial  carbon  temperature  =  740  F. 

4.   CASE  V-4   Reaction  Order  =  1.00 

For  this  case,  the  temperature  of  the  combustion  and 
extinction  processes  are: 


for  combustion  Tc  =  670  F 


for  extinction  T. 


660  F 


The  Figures  5.19  to  5.24  show  the  behavior  of  the 
carbon  temperature,  oxygen  concentration  and  reaction  at 
each  position  and  time. 

D .   SUMMARY 

1 .   Power  Relation 

For  each  case,  a  pair  of  values  (reaction  order, 
combustion  temperature)  was  obtained. The  points  were  plotted 
in  cartesian  coordinates  in  Figure  5.25  .  A  log-log  plot  is 
shown  in  Figure  5.26  ,  and  a  semi  log  plot  is  shown  in 
Figure  5.27  .  It  is  observed  that  the  relation  between 
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TEMPERATURE  SURFACE  FROM  GRAF3E 


Figure  5.19    Temperature  vs  X/L  and  time 

for  reaction  order  =  1.00 

Initial  carbon  temperature  =  670  F. 


reaction  order  (n)  and  combustion  temperature  (T  )  of  figure 
5.27  yields  the  approximate  power  relation  (equation  5.2)  . 
Table  5.25  shows  a  resume  of  each  case,  combustion 
temperature,  extinction  temperature  and  the  value  of  the 
combustion  temperature  obtained  from  Equation  5.2  and  the 
percentage  error  .  An  observation  of  these  results  leads  to 
the  conclusion  that  this  power  relation  is  a  good  estimation 
for  combustion  temperature. 
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temperature:  surface  from  graf3E 


Figure  5.20    Temperature  vs  X/L  and  time 

for  reaction  order  =  1.00 

Init-ial  carbon  temperature  =  660  F. 

OXYGEN  CONC.  SURFACE  FROM  GRAF3E 
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Figure  5.21    Oxygen  concentration  vs  X/L  and  time 

for  reaction  order  =  1.00 

Initial  carbon  temperature  =  670  F. 
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OXYGEN  CONC.  SURFACE  FRCM  GRAr3E 


Figure  5.22    Oxygen  concentration  vs  X/L  and  time 

for  reaction  order  =  1.00 

Initial  carbon  temperature  =  660  F. 
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Figure  5.23    Reaction  rate  vs  X/L  and  time 

for  reaction  rate=  1.00 

Initial  carbon  temperature  =  670  F. 
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RERCTION  RATE  SURFRCE  FROM  GRRF3E 


T"«  nee. 


Figure  5.24    Reaction  rate  vs  X/L  and  time 

for  reaction  order  =  1.00 

Initial  carbon  temperature  =  660  F. 


T   =  L0(305-022n)   (  F) 


(5.2) 
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reaction 

program 

order 

output 

0.25 

970 

0.50 

840 

0.75 

740 

1.00 

670 

TABLE  XV 

COMBUSTION  AND  EXTINCTION  TEMPERATURE 
FOR  REACTION  ORDER  CASES 

temperatures  (QF) 
extinction      combustion  %  error 


0.0 

1.6 
1.4 
0.0 


program 

power 

output 

relation 

980 

980 

850 

863 

750 

760 

670 

670 
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REACTION  ORDER  VERSUS  TEMPERATURE 
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Figure  5.25    Reaction  order  versus  temperature 

Retangular  plot. 
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REACTION  ORDER  VERSUS  TEMPERATURE 


2M0" 


REACTION  ORDER 


Figure  5.26    Reaction  order  versus  temperature 

Log- log  plot . 
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REACTION  ORDER  VERSUS  TEMPERATURE 
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Figure  5.27 


Reaction  order  versus  temperature 
Semi  log  plot. 
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VI.  CONCLUSIONS 

In  each  of  the  previous  chapters,  the  effects  of  a 
single  parameter  on  the  combustion  behavior  of  a  porous 
medium  was  investigated.  In  Chapter  2,  we  considered  the 
effect  of  heat  flux  (SQ)  and  its  duration  (TQ)  on 
combustion.  In  Chapter  3,  the  effect  of  medium  thickness  (L) 
on  behavior  was  investigated.  Chapter  4  was  a  study  of  the 
effect  of  permeability  (m)  on  system  behavior;  and  finally 
in  Chapter  5,  the  effect  of  reaction  order  (n)  was 
investigated . 

The  major  conclusions  from  these  investigations  are: 

-  The  relation  between  SQ  and  TQ  is  not  linear  ,  a  power 
relation  gives  rasonable  approximation. 

-  During  the  combustion  process,  for  the  SQ  cases, 
initially  there  is  a  large  temperature  gradient  across 
the  medium  with  T  (  X/L  =  0 . )  »  T  (  X/L  =  1.)  .  With 
passing  time,  this  gradient  decreases  until  there  is  a 
uniform  combustion  temperature  over  the  entire  medium. 
When  this  occurs ,  there  is  no  oxygen  within  the 
medium,  and  the  medium  is  continuosly  consumed  by 
combustion . 

-  The  relation  between  thickness  and  combustion  tempera 
ture  is  not  linear  .  Increasing  the  thickness  results 
in  decreasing  of  combustion  temperature.  If  the 
thickness  increased  with  a  constant  rate  c  (  Lj+i/1j  = 
c  )  the  time  that  marks  the  start  of  combustion  is  not 
constant  ,  it  increases  with  increasing  thickness  ( 
i.e.  ti+i/t\  "  t»(L)  )  where  b  depends  upon  the 
thickness  . 

-  The  relation  between  permeability  and  minimum  combus 
tion  temperature  is  not  linear. 

-  An  increase  of  the  permeability  shows  that  the  combus 
tion  temperature  increases  to  an  asymptotic  limit. 
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-  Permeability  has  a  significant  influence  on  the  speed 
of  the  combustion  process.  When  the  permeability  is 
increased,  the  time  that  marks  the  start  of  the 
combustion  process  decreases.  For  large  values  of 
permeability,  the  •  combustion  may  occur  in  a  fraction 
of  a  second. 

-  If  the  permeability  is  increased  with  a  constant  rate 
c(  nij+-i  =  c  m^  )  ,  the  time  that  measures  the  start  of 
combustion  decreases  . 

-  The  relation  between  reaction  order  (n)  and  combustion 
temperature  (Tc)  is  not  linear. 

-  An  observation  of  the  graphical  surfaces,  for  all 
cases  of  this  investigation  ,  shows  a  similar  form  for 
all  combustion  and  extinction  cases.  For  all 
combustion  cases,  it  was  noted  that  the  oxygen  within 
the  porous  medium  being  consumed  with  passing  time. 
The  first  point  at  which  the  oxygen  is  depleted  is 
always  X/L  =  1.  With  passing  time  the  point  with 
smaller  X/L  reach  zero  oxygen,  until  all  points  have 
zero  oxygen.  At  this  time  the  only  oxygen  is  at  the 
interface  of  the  porous  medium  and  the  ambient 
environment.  At  this  time  the  combustion  process 
becomes  a  surface  recession  phenomena.  At  the  same 
time  the  oxygen  is  being  depleted  the  temperature  is 
increasing.  Eventually  when  all  the  interior  oxygen 
has  been  consumed,  there  is  a  temperature  gradient 
across  the  medium  with  Tc(X/L=0)  >  T  (X/L=l). 

This  thesis  concludes  with  two  recomendations  .  First  , 
it  is  recommended  that  additional  parameter  studies  he 
conducted  .  For  example  ,  it  would  be  useful  to  determine 
the  effects  of  pressure  differential  ,  conductivity  ,  and 
other  thermophysical  properties  on  combustion  .  Another 
recommendation  is  to  attempt  to  obtain  a  relation  between 
combustion  temperature  and  several  parameters  . 
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